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Abstract 

These lectures contain an introduction to supersymmetric theories and the minimal su- 
persymmetric standard model. Phenomenological and cosmological consequences of super- 
symmetry are also discussed. 



1 Based on lectures delivered at the Les Houches Summer School, July 1999 



1 Introduction 



It is often called the last great symmetry of nature. Rarely has so much effort, both theoret- 
ical and experimental, been spent to understand and discover a symmetry of nature, which 
up to the present time lacks concrete evidence. Hopefully, in these lectures, where I will give 
a pedagogical description of supersymmetric theories, it will become clear why there is so 
much excitement concerning supersymmetry's role in nature. 

After some preliminary background on the standard electroweak model, and some mo- 
tivation for supersymmetry, I will introduce the notion of supersymmetric charges and the 
supersymmetric transformation laws. The second lecture will present the simplest supersym- 
metric model (the non-interacting massless Wess-Zumino model) and develop the properties 
of chiral superfields, auxiliary fields, the superpotential, gauge multiplets and interactions. 
The next two lectures focus on the minimal supersymmetric standard model (MSSM) and 
its constrained version which is motivated by supergravity. The last two lectures will look 
primarily at the cosmological and phenomenological consequences of supersymmetry. 

1.1 Some Preliminaries 

Why Supersymmetry? If for no other reason, it would be nice to understand the origin of 
the fundamental difference between the two classes of particles distinguished by their spin, 
fermions and bosons. If such a symmetry exists, one might expect that it is represented by 
an operator which relates the two classes of particles. For example, 



As such, one could claim a raison d'etre for fundamental scalars in nature. Aside from the 
Higgs boson (which remains to be discovered), there are no fundamental scalars known to 
exist. A symmetry as potentially powerful as that in eq. (P is reason enough for its study. 
However, without a connection to experiment, supersymmetry would remain a mathemat- 
ical curiosity and a subject of a very theoretical nature as indeed it stood from its initial 
description in the early 1970's fl], |2| until its incorporation into a realistic theory of physics 
at the electroweak scale. 

One of the first break-throughs came with the realization that supersymmetry could help 
resolve the difficult problem of mass hierarchies ||, namely the stability of the electroweak 
scale with respect to radiative corrections. With precision experiments at the electroweak 
scale, it has also become apparent that grand unification is not possible in the absence of 
supersymmetry [|J . These issues will be discussed in more detail below. 

Because one of our main goals is to discuss the MSSM, it will be useful to first describe 
some key features of the standard model if for no other reason than to establish the notation 
used below. The standard model is described by the SU(3) c x SU(2) L x U(l)y gauge group. 
For the most part, however, I will restrict the present discussion to the electroweak sector. 
The Lagrangian for the gauge sector of the theory can be written as 



Q | Boson) 
Q | Fermion) 



| Fermion) 
| Boson) 



(1) 



(2) 
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where = d^W* — d u W^ + ge^ k W 3 J¥ k is the field strength for the SU(2) gauge boson 
WL and F^ u = d^B u — d u B^ is the field strength for the U{1) gauge boson £> M . The fermion 
kinetic terms are included in 

^F = ~J2 i [fL-fDnfL + ftffDfJn] (3) 

/ 

where the gauge covariant derivative is given by 



D^dp-ig^Wi-igf^Bp (4) 



The (7j are the Pauli matrices (representations of SU(2)) and Y is the hypercharge. g and g' 
are the SU(2) L and U(l)y gauge couplings respectively 

Fermion mass terms are generated through the coupling of the left- and right-handed 
fermions to a scalar doublet Higgs boson 0. 

CY = -Y,[Gf<j>fLf R ]+h.c. (5) 
/ 

The Lagrangian for the Higgs field is 

£^ = -|ZV>| 2 -V(0) (6) 
where the (unknown) Higgs potential is commonly written as 

V(4>) = V0V + A(0W (7) 
The vacuum state corresponds to a Higgs expectation valuef] 

(0> = <0*>=(°) with v 2 = ^ (8) 

The non-zero expectation value and hence the spontaneous breakdown of the electroweak 
gauge symmetry generates masses for the gauge bosons (through the Higgs kinetic term in 
(^) and fermions (through (|5[)). In a mass eigenstate basis, the charged H^-bosons (W = 
(W 1 ± iW 2 )/\/2) receive masses 

1 



M w = -^gv (9) 



The neutral gauge bosons are defined by 



_ gWj - g'B, g'Wj + gB, 



with masses 

1 



M z = -=^g 2 + g ,2 v = M w /cose w m 7 = (11) 



2 Note that the convention used here differs by a factor of \[2 from that in much of the standard model 
literature. This is done so as to conform with the MSSM conventions used below. 
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where the weak mixing angle is defined by 



sme w = g'/yjg 2 + g' 2 (12) 

Fermion masses are 

rrif = GfV (13) 

As one can see, there is a direct relationship between particle masses and the Higgs 
expectation value, v. Indeed, we know from (^) and ([□]) that v ~ M w ~ 0(100) GeV. We 
can then pose the question, why is Mw <C Mp = 1.2 x 10 19 GeV or equivalently why is 



1.2 The hierarchy problem 

The mass hierarchy problem stems from the fact that masses, in particular scalar masses, are 
not stable to radiative corrections ||. While fermion masses also receive radiative corrections 
from diagrams of the form in Figure [l], these are only logarithmically divergent (see for 
example H), 

5rrif ~ — nif ln(A 2 /m 2 ) (14) 

A is an ultraviolet cutoff, where we expect new physics to play an important role. As one 
can see, even for A ~ M P , these corrections are small, 5mj < m,f. 




Figure 1: 1-loop correction to a fermion mass. 

In contrast, scalar masses are quadratically divergent. 1-loop contributions to scalars 
masses, such as those shown in Figure || are readily computed 

5m 2 H ~ g}, g\ X J d 4 k± ~ 0(£)A 2 (15) 

due to contributions from fermion loops with coupling gf, from gauge boson loops with 
coupling g 2 , and from quartic scalar-couplings A. From the relation @ and the fact that the 
Higgs mass is related to the expectation value, m 2 H = 4t> 2 A, we expect Mw ~ ran- However, 
if new physics enters in at the GUT or Planck scale so that A ^> Mw, the 1-loop corrections 
destroy the stability of the weak scale. That is, 

A > M w -> 5m 2 H > m 2 H (16) 

Of course, one can tune the bare mass m# so that it contains a large negative term 
which almost exactly cancels the 1-loop correction leaving a small electroweak scale mass 2 . 
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For a Planck scale correction, this cancellation must be accurate to 32 significant digits. 
Even so, the 2-loop corrections should be of order a 2 A 2 so these too must be accurately 
canceled. Although such a series of cancellations is technically feasible, there is hardly a 
sense of satisfaction that the hierarchy problem is under control. 

An alternative and by far simpler solution to this problem exists if one postulates that 
there are new particles with similar masses and equal couplings to those responsible for the 
radiatively induced masses but with a difference (by a half unit) in spin. Then, because 
the contribution to 5m 2 H due to a fermion loop comes with a relative minus sign, the total 
contribution to the 1-loop corrected mass 2 is 

6m% ~ 0(^)(A 2 + ml) - 0(£)(A 2 + m F ) = 0(£)(m| - m 2 F ) (17) 

If in addition, the bosons and fermions all have the same masses, then the radiative correc- 
tions vanish identically. The stability of the hierarchy only requires that the weak scale is 
preserved so that we need only require that 

\m 2 B -m 2 F \ < 1 TeV 2 (18) 

As we will see in the lectures that follow, supersymmetry offers just the framework for 
including the necessary new particles and the absence of these dangerous radiative corrections 

Before we embark, I would like to call attention to some excellent additional resources on 
supersymmetry. These are the classic by Bagger and Wess on supersymmetry, [0, the book 



by Ross on Grand Unification || and two recent reviews by Martin || and Ellis |T0] 



1.3 Supersymmetric operators and transformations 

Prior to the introduction of supersymmetry, operators were generally regarded as bosonic. 
That is, they were either scalar, vector, or tensor operators. The momentum operator, P M , is 
a common example of a vector operator. However, the types of bosonic charges are greatly 



limited, as was shown by Coleman and Mandula Given a tensorial operator, E^, its 
diagonal matrix elements can be decomposed as 

< a|E Mi/ |a >= atpffi + PQuv (19) 
One can easily see that unless a = 0, 2 to 2 scattering process allow only forward scattering. 
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Operators of the form expressed in (ffl) however, are necessarily non-diagonal as they 
require a change between the initial and final state by at least a half unit of spin. Indeed, 
such operators, if they exist must be fermionic in nature and carry a spinor index Q a . There 
may in fact be several such operators, Q l a with i — 1, ...,N, (though for the most part we 
will restrict our attention to N — 1 here). As a symmetry operator, Q must commute with 
the Hamiltonian H, as must its anti-commutator. So we have 

[Q^H} = [{<£,<#}, #] = (20) 

By extending the Coleman- Mandula theorem ||12|| , one can show that 

{Q\Q j] } oc fiPp + ZF (21) 

where Z l i is antisymmetric in the supersymmetry indices {i, j}. Thus, this so-called "central 
charge" vanishes for N = 1. More precisely, we have in a Weyl basis 

{Q Q ,Qy = 2a^P, 

{Q a ,Qp} = {Ql,Q^} = 

[Q ai P,] = [Q^,P,] = (22) 

Before setting up the formalism of the supersymmetric transformations, it will be useful to 
establish some notation for dealing with spinors in the Dirac and Weyl bases. The Lagrangian 
for a four- component Dirac fermion with mass M, can be written as 



where 



C D = -Mo^d^D - M^d^d (23) 

*) -to 1 \) M 

and = (l,(Xj), a M = (1, — <7j), Oi are the ordinary 2x2 Pauli matrices. I am taking the 
Minkowski signature to be (—,+,+,-1-). We can write the Dirac spinor ^> £, in terms of 2 
two-component Weyl spinors 

*D = U a e f «) (25) 

Note that the spinor indices (a, a) are raised and lowered by where {ij} can be either 
both dotted or both undotted indices, e is totally antisymmetric and = — e* J with e 12 = 1. 
It is also useful to define projection operators, Pl and Pr with 

with a similar expression for Pr. In this way we can interpret £ a as a left-handed Weyl spinor 
and as a right-handed Weyl spinor. The Dirac Lagrangian (p3|) can now be written in 
terms of the two-component Weyl spinors as 

Cd = -i^d^ - ix^d^x ~ M{i X + £V) (27) 
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having used the identity, — x alx £, = C^X- 

Instead, it is sometimes convenient to consider a four-component Majorana spinor. This 
can be done rather easily from the above conventions and taking £ = x, so that 

*"=(J*) Vm = (C (28) 

and the Lagrangian can be written as 

i— 1 — 

= -it}o»d^- ^M(ee+et£t) (29 ) 

The massless representations for supersymmetry are now easily constructed. Let us 
consider here N — 1 supersymmetry, i.e., a single supercharge Q a . For the massless case, 
we can choose the momentum to be of the form P M = |(— 1, 0, 0, 1). As can be readily found 
from the anticommutation relations fl2"2|), the only non- vanishing anticommutation relation 
is {Qi,Q\} = 1. Consider then a state of given spin, | A) such that Q\\\) = 0. (If it is not 
0, then due to the anticommutation relations, acting on it again with Q\ will vanish.) From 
the state |A), it is possible to construct only one other nonvanishing state, namely Qi\\) - 
the action of any of the other components of Q a will vanish as well. Thus, if the state |A) 
is a scalar, then the state Qi|A) will be a fermion of spin 1/2. This (super)multiplet will be 
called a chiral multiplet. If |A) is spin 1/2, then Qi|A) is a vector of spin 1, and we have a 
vector multiplet. In the absence of gravity (supergravity), these are the only two types of 
multiplets of interest. 

For N > 1, one can proceed in an analogous way. For example, with N = 2, we begin 
with two supercharges Q l ,Q 2 . Starting with a state |A), we can now obtain the following: 
Qi|A), Q\\\)i <5iQi|A). In this case, starting with a complex scalar, one obtains two fermion 
states, and one vector, hence the vector (or gauge) multiplet. One could also start with a 
fermion state (say left-handed) and obtain two complex scalars, and a right-handed fermion. 
This matter multiplet however, is clearly not chiral and is not suitable for phenomenology. 
This problem persists for all supersymmetric theories with N > 1, hence the predominant 
interest in N = 1 supersymmetry. 

Before we go too much further, it will be useful to make a brief connection with the 
standard model. We can write all of the standard model fermions in a two-component Weyl 
basis. The standard model fermions are therefore 
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Li 



'Li 



t*L, T L (30) 



Note that the fields above are all left-handed. Color indices have been suppressed. From 
(p9j). we see that we would write the fermion kinetic terms as 

Ckin = -iQla^Qi - mfWdpuZ (31) 



* e =(S) = f! L ) (32) 



As indicated above and taking the electron as an example, we can form a Dirac spinor 

J L / \ e R 

A typical Dirac mass term now becomes 

^ e ^ e = e c L e L + e\e c j} = e ] R e L + e[e R (33) 

As we introduce supersymmetry, the field content of the standard model will necessarily be 
extended. All of the standard model matter fields listed above become members of chiral 
multiplets in N = 1 supersymmetry. Thus, to each of the (Weyl) spinors, we assign a 
complex scalar superpartner. This will be described in more detail when we consider the 
MSSM. 

To introduce the notion of a supersymmetric transformation, let us consider an infinites- 
imal spinor £ a with the properties that £ anticommutes with itself and the supercharge Q, 
but commutes with the momentum operator 

{c , = {e , Qp) = [p» r] = o (34) 

It then follows that since both £ and Q commute with P M , the combination £Q also commutes 
with or 

[P l *,ZQ] = [P fl ,Zty] = (35) 

where by £Q we mean £Q = £ a Q a = e a ^ a Q p = -e af3 Q^ a = e^ a Q^ a = Qf . Similarly, 

= ^aQ^ ■ Also note that (, a Q a = —£, a Q a - Finally, we can compute the commutator of 
£gand£tQt ) 

[£Q, ?Q*] = iQi ] Q ] - £ ] Q ] £Q = CQa^&Q^ - i\Q ]P CQ a 



= zfr^Pv (36) 

We next consider the transformation property of a scalar field, 0, under the infinitesimal 

W = (CQa + ^ V (37) 
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As described above, we can pick a basis so that Q^(j) = 0. Let call the spin 1/2 fermion 

Q a (p = \/2V>a, so 

Srf = CQa<P = V2C^ a (38) 

To further specify the supersymmetry transformation, we must define the action of Q and 
on ip. Again, viewing Q as a "raising" operator, we will write Qa^-y = ~V2e ai F and 

Q^ip-y = —y/2ia^dfj,(j), where F, as we will see, is an auxiliary field to be determined later. 
Even though we had earlier argued that Q acting on the spin 1/2 component of the chiral 
multiplet should vanish, we must keep it here, though as we will see, it does not correspond 
to a physical degree of freedom. To understand the action of Q\ we know that it must be 
related to the scalar <ft, and on dimensional grounds (Q^X is of mass dimension 3/2) it must 
be proportional to P^cj). Then 

= (r q q + 6Q f V 7 

= -V2e £oer F + V2i^%((> 

= v^^ + ^K^)^ ( 39 ) 

Given these definitions, let consider the successive action of two supersymmetry trans- 
formations 5% and 8 V . 

= V2[-V2C a v% a F + V2iC a a^%(f>} (40) 

If we take the difference (S V S^ — 5^5 r) )(p, we see that the first term in ( fiO|) cancels if we write 
i a rpe^ a = —£, a rj a and note that (, a r] a = r] a ^ a . Therefore the difference can be written as 

(V? - = 2fa<7"£ t - £<rV W ( 41 ) 

In fact it is not difficult to show that fl41~D is a result of the general operator relation 

6 v 6t - 5s5 v = 2( W ^ - (42) 

Knowing the general relation ( |4"2]) and applying it to a fermion r/> 7 will allow us to deter- 
mine the transformation properties of the auxiliary field F. Starting with 

S v 5^ = -V2Ce ai 5 v F + 2ia^ V a d^ a (43) 

we use the Fierz identity X-y^Ca + V-y( a Xa + C-yX^a = 0, and making the substitutions, 
Xi = rfp&^i V = V, and ( = d^, we have 

5 V S^ = -V2Ce ai 6 v F - 2i{r] 7 d^ a a^ + d^a»% ^ Va } (44) 

Next we use the spinor identity, x a £^ — — & p ~o^ a Xa along with rpXi — X 7? ?7 (from 
above) to get 

VfV>7 = -V2C a e ai 5 v F + 2i{r h ? $ a^ a d^ a - V a a^d^} (45) 
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It is not hard to see then that the difference of the double transformation becomes 
(Vs-W^ = 2(r^ t -W)^ 7 

+V2{^6 V F - rjyS^F) (46) 
Thus, the operator relation ( fl^) will be satisfied only if 

5^F = -V2(^P^) (47) 
and we have the complete set of transformation rules for the chiral multiplet. 

2 The Simplest Models 

2.1 The massless non-interacting Wess-Zumino model 

We begin by writing down the Lagrangian for a single chiral multiplet containing a complex 
scalar field and a Majorana fermion 

= -d^*dy- % ~{^d^-d tl ^a^) (48) 

where the second line in ([4^) is obtained by a partial integration and is done to simplify the 
algebra below. 

We must first check the invariance of the Lagrangian under the supersymmetry transfor- 
mations discussed in the previous section. 

+-j=£%F*cr> t il> + ^a^d^cj)* + h.c. (49) 

Now with the help of still one more identity, (a fJ- a 1 ' + o^cx^)^ = — 2r] tiu 5^ , we can expand the 
above expression 

-J-g^o^d^ - d„F*a^) + h.c. (50) 
v2 

Fortunately, we now have some cancellations. The first and third terms in (|50|) trivially 
cancel. Using the commutivity of the partial derivative and performing a simple integration 
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by parts we see that the second and fourth terms also cancel. We left with (again after an 
integration by parts) 

5C = -iV^F^d^ + h.c. (51) 
indicating the lack of invariance of the Lagrangian (ESI). 



We can recover the invariance under the supersymmetry transformations by considering 
in addition to the Lagrangian (f48|) the following, 



C aux = F*F (52) 

and its variation 

SC aux = SF*F + F*SF (53) 



The variation of the auxiliary field, F, was determined in fl47|) and gives 

5C aux = iV2gF*a»d^ + h.c. (54) 
and exactly cancels the piece left over in fl5"l~|). Therefore the Lagrangian 

C = -d^d^ - i^a^ip + F*F (55) 
is fully invariant under the set of supersymmetry transformations. 

2.2 Interactions for Chiral Multiplets 

Our next task is to include interactions for chiral multiplets which are also consistent with 
supersymmetry. We will therefore consider a set of chiral multiplets, ((pi, ipi, Fi) and a renor- 
malizable Lagrangian, Ci nt . Renormalizability limits the mass dimension of any term in the 
Lagrangian to be less than or equal to 4. Since the interaction Lagrangian must also be 
invariant under the supersymmetry transformations, we do not expect any terms which are 
cubic or quartic in the scalar fields 0j. Clearly no term can be linear in the fermion fields 
either. This leaves us with only the following possibilities 

C int = -zAHtbi + B l Fi + h.c. (56) 

where is some linear function of the <pi and (f> 1 * and B l is some function which is at most 
quadratic in the scalars and their conjugates. Here, and in all that follows, it will be assumed 
that repeated indices such as ii are summed. Furthermore, since ipiipj = ipjipi (spinor indices 
are suppressed), the function A % i must be symmetric in ij. As we will see, the functions A 
and B will be related by insisting on the invariance of (|56|) . 
We begin therefore with the variation of Ci n t 

IdAv , /- , 1 dA ij 



I / 1 /t " j I / 1 /i "j i 
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+B i V2i^a' x d fl iP i + h.c. (57) 

where the supersymmetry transformations of the previous section have already been per- 
formed. The notation (ipiipj) refers to ipfipj as clearly spinor indices have everywhere been 
suppressed. The Fierz identity (£ipk)('4 ) i' l l>j) + (£V'i)('^jV'fc) + (£'0.j)(V'fc'0i) — implies that 
the derivative of the function A l i with respect to <pk (as in the first term of Q57D ) must be 
symmetric in ijk. Because there is no such identity for the second term with derivative with 
respect to <p k , this term must vanish. Therefore, the function A % i is a holomorphic function 
of the 4>i only. Given these constraints, we can write 

A ij = _ M H _ y ijk ( j )k ( 58 ) 

where by (|56|) we interpret M y as a symmetric fermion mass matrix, and y y as a set of 
(symmetric) Yukawa couplings. In fact, it will be convenient to write 

d 2 W 

where 

w = lM ij <f> i( f>i + ly lJk MAk (60) 

2 6 

and is called the superpotential. 

Noting that the 2nd and 3rd lines of (|57|) are equal due to the symmetry of we can 
rewrite the remaining terms as 

5C mt = A^^^Fj + V^ia^dM 

+ _y 2mFi + w( ^ )Fi 

-B^id^a^ + h.c. (61) 

using in addition one of our previous spinor identities on the last line. Further noting that 
because of our definition of the superpotential in terms of A % \ we can write A^d^j = 
—d^dW/dtpi). Then the 2nd and last terms of (|61~1) can be combined as a total derivative if 

dW 

d(j)i 

and thus is also related to the superpotential W . Then the 4th term proportional to dB l / dft 
is absent due to the holomorphic property of XV, and the definition of B (|62]) allows for a 
trivial cancellation of the 1st and 3rd terms in Q6T|). Thus our interaction Lagrangian 
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is in fact supersymmetric with the imposed relationships between the functions A l i , B l , and 
the superpotential W. 

After all of this, what is the auxiliary field Fl It has been designated as an "auxiliary" 
field, because it has no proper kinetic term in (|55|) . It can therefore be removed via the 
equations of motion. Combining the Lagrangians (|55| ) and ([51]) we see that the variation of 
the Lagrangian with respect to F is 

5 A = F l * + W i (63) 
or 

where we can now use the convenient notation that W l = dW/d(pi, W* = dW/dcp 1 *, and 
yyij _ Q 2 W I Q^fi^).^ e tc. The vanishing of (|63"D then implies that 

Fi = -W* (64) 

Putting everything together we have 

+ W*^^ j ) - W l W* (65) 

As one can see the last term plays the role of the scalar potential 

V{(j) h ^) = W i W* (66) 



2.3 Gauge Multiplets 

In principle, we should next determine the invariance of the Lagrangian including a vector 
or gauge multiplet. To repeat the above exercise performed for chiral multiplets, while 
necessary, is very tedious. Here, I will only list some of the more important ingredients. 

Similar to the definition in (£|), the gauge covariant derivative acting on scalars and chiral 
fermions is 

Dp — — igT ■ (67) 

where T is the relevant representation of the gauge group. For SU(2), we have simply that 
T % = cr*/2. In the gaugino kinetic term, the covariant derivative becomes 

(IV) a = d M A a + £/^7l 6 A c (68) 

where the f abc are the (antisymmetric) structure constants of the gauge group under con- 
sideration ([T a ,T b ] = if abc T c ). Gauge invariance for a vector field, A a ^ is manifest through 
a gauge transformation of the form 

5„A; = -d,A a + gf^A'Al (69) 

where A is an infinitesimal gauge transformation parameter. To this, we must add the gauge 
transformation of the spin 1/2 fermion partner, or gaugino, in the vector multiplet 

S„\ a = gf abc A b \ c (70) 
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Given our experience with chiral multiplets, it is not too difficult to construct the super- 
symmetry transformations for the the vector multiplet. Starting with A^, which is real, and 
taking QA^ = 0, one finds that 

6 e Al = i[£a»\ a - \ a %£] (71) 

Similarly, applying the supersymmetry transformation to A a , leads to a derivative of A^ (to 
account for the mass dimension) which must be in the form of the field strength F®, and 
an auxiliary field, which is conventionally named D a . Thus, 

5^ a = l -{a^i)F; u + iiD a (72) 

As before, we can determine the transformation properties of D a by applying successive 

supersymmetry transformations as in ( pE2"D with the substitution <9 M — > using Q57| ) above. 
The result is, 

5^D a = [ga> i D l X + £> M A°Vf ] (73) 

Also in analogy with the chiral model, the simplest Lagrangian for the vector multiplet 

is 

Cgauge = ~F^F^ - iX^D^ + JlW (74) 



In (|74]), the gauge kinetic terms are given in general by = A" — d u A® + gf abc A b ^A c u . 

If we have both chiral and gauge multiplets in the theory (as we must) then we must 
make simple modifications to the supersymmetry transformations discussed in the previous 
section and add new gauge invariant interactions between the chiral and gauge multiplets 
which also respect supersymmetry. To (p9|), we must only change <9 M — > using fl6"T|). To 
([47l), it will be necessary to add a term proportional to {T a (j))^X a ^ so that, 

5^F = v^V/^) + 2^(T»£ t A at (75) 

The new interaction terms take the form 

C mt = V2gi [(0*T»A a - At a (^T»l + g{<P*T a <P)D a (76) 



Furthermore, invariance under supersymmetry requires not only the additional term in (|75|) , 
but also the condition 

\V : (T')',o :i = (77) 

Finally, we must eliminate the auxiliary field D a using the equations of motion which 
yield 

D a = -g{<t>*T a (j>) (78) 
Thus the "D-term" is also seen to be a part of the scalar potential which in full is now, 

V(<f), 0*) = \Ff + ^\D a \ 2 = \Wf + i/(0*T» 2 . (79) 

Notice a very important property of the scalar potential in supersymmetric theories: the 
potential is positive semi-definite, V > 0. 
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2.4 Interactions 

The types of gauge invariant interactions allowed by supersymmetry are scattered throughout 
the pieces of the supersymmetric Lagrangian. Here, we will simply identify the origin of a 
particular interaction term, its coupling, and associated Feynmann diagram. In the diagrams 
below, the arrows denote the flow of chirality. Here, tjj will represent an incoming fermion, 
and ijfl an outgoing one. While there is no true chirality associated with the scalars, we 
can still make the association as the scalars are partnered with the fermions in a given 
supersymmetric multiplet. We will indicate <fi with an incoming scalar state, and <p* with an 
outgoing one. 

Starting from the superpotential flB0| ) and the Lagrangian we can identify several 
interaction terms and their associated diagrams: 

• a fermion mass insertion from W^ip^j 



• a scalar-fermion Yukawa coupling, also from W^i^^j 



k 




• a scalar mass insertion from IW 1 ] 2 
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a scalar cubic interaction from |W^| 2 (plus its complex conjugate which is not shown) 







s J 



< 



s k 



• and finally a scalar quartic interaction from \ W 



i\2 



\ / 
\ / 



k 



\ 



\ 

j ' i 



Next we must write down the interactions and associated diagrams for the gauge mul- 
tiplets. The first two are standard gauge field interaction terms in any non-abelian gauge 
theory (so that f abc 7^ 0) and arise from the gauge kinetic term F^ u , the third is an interac- 
tion between the vector and its fermionic partner, and arises from the gaugino kinetic term 
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• The quartic gauge interaction from F^ u (to be summed over the repeated gauge indices) 



,2 jabc jade j^bV j^cv j^d j^e 




• The trilinear gauge interaction also from 



gfobc A b^ A cu^ A a v _ Q^aj 




• The gauge-gaugino interaction from X^a^D^X 



gfo-bc A b^ X a^ X c 




If our chiral multiplets are not gauge singlets, then we also have interaction terms between 
the vectors and the fermions and scalars of the chiral multiplet arising from the chiral kinetic 
terms. Recalling that the kinetic terms for the chiral multiplets must be expressed in terms 
of the gauge covariant derivative fl6"7|), we find the following interactions, from|D M 0| 2 and 
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• a quartic interaction interaction involving two gauge bosons and two scalars, 



g 2 A afl A b ll (T a ( f))((j)*T b ) 




• a cubic interaction involving one gauge boson and two scalars, 



giA^iT^d^* + h.c.) vAAAAA 



a cubic interaction involving one gauge boson and two fermions, 



g(A a ^a»(T a il>) + h.c.) 




Finally, there will be two additional diagrams. One coming from the interaction term 
involving both the chiral and gauge multiplet, and one coming from D 2 , 
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• a cubic interaction involving a gaugino, and a chiral scalar and fermion pair, 



g{{(t)*T a i))\ a + h.c.) 



another quartic interaction interaction involving a gaugino, and a chiral scalar and 
fermion pair, 

\ / 
\ / 

\ / 

/ \ 



2.5 Supersymmetry Breaking 

The world, as we know it, is clearly not supersymmetric. Without much ado, it is clear from 
the diagrams above, that for every chiral fermion of mass M, we expect to have a scalar 
superpartner of equal mass. This is, however, not the case, and as such we must be able to 
incorporate some degree of supersymmetry breaking into the theory. At the same time, we 
would like to maintain the positive benefits of supersymmetry such as the resolution of the 
hierarchy problem. 

To begin, we must be able to quantify what we mean by supersymmetry breaking. From 
the anti-commutation relations (^), we see that we can write an expression for the Hamil- 
tonian or P° using the explicit forms of the Pauli matrices as 

^ = 7 <&} (80) 

A supersymmetric vacuum must be invariant under the supersymmetry transformations and 
therefore would require Q\0) = and Q^\0) = and therefore corresponds to H = and 
also V = 0. Thus, the supersymmetric vacuum must have \F\ = \D\ = 0. Conversely, if su- 
persymmetry is spontaneously broken, the vacuum is not annihilated by the supersymmetry 
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charge Q so that Q\0) = x an d (x|Q|0) — where \ is a fermionic field associated with the 
breakdown of supersymmetry and in analogy with the breakdown of a global symmetry, is 
called the Goldstino. For f x ^ 0, (0\H\0) — Vo ^ 0, and requires therefore either (or both) 
\F\ 7^ or \D\ 7^ 0. Mechanisms for the spontaneous breaking of supersymmetry will be 
discussed in the next lecture. 

It is also possible that to a certain degree, supersymmetry is explicitly broken in the 
Lagrangian. In order to preserve the hierarchy between the electroweak and GUT or Planck 
scales, it is necessary that the explicit breaking of supersymmetry be done softly, i.e., by the 
insertion of weak scale mass terms in the Lagrangian. This ensures that the theory remain 
free of quadratic divergences The possible forms for such terms are 



where the are gaugino masses, m 2 are soft scalar masses, B is a bilinear mass term, and 
A is a trilinear mass term. Masses for the gauge bosons are of course forbidden by gauge 
invariance and masses for chiral fermions are redundant as such terms are explicitly present 
in M !J already. The diagrams corresponding to these terms are 

• a soft supersymmetry breaking gaugino mass insertion 




l -{BMfMj ~ ^{Ayf^^u + h.c. 



(81) 



M?A a A 



x 



x 



a 



-X 



• a soft supersymmetry breaking scalar mass insertion 




• a soft supersymmetry breaking bilinear mass insertion 
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a soft supersymmetry breaking trilinear scalar interaction 



y J 



We are now finally in a position to put all of these pieces together and discuss realistic 
super symmetric models. 

3 The Minimal Supersymmetric Standard Model 



To construct the supersymmetric standard model [14| we start with the complete set of chiral 



fermions in (|30|) , and add a scalar superpartner to each Weyl fermion so that each of the fields 
in fl30|) represents a chiral multiplet. Similarly we must add a gaugino for each of the gauge 
bosons in the standard model making up the gauge multiplets. The minimal supersymmetric 
standard model (MSSM) ]T3] is defined by its minimal field content (which accounts for the 
known standard model fields) and minimal superpotential necessary to account for the known 
Yukawa mass terms. As such we define the MSSM by the superpotential 

W = eij [y e H{Ve c + y d H{Q i d c + y u H^u c ] + W» (82) 

where 

W, = e^HlHi (83) 

In (^), the indices, {ij}, are SXJ(2) L doublet indices. The Yukawa couplings, y, are all 3 x 3 
matrices in generation space. Note that there is no generation index for the Higgs multiplets. 
Color and generation indices have been suppressed in the above expression. There are two 
Higgs doublets in the MSSM. This is a necessary addition to the standard model which can 
be seen as arising from the holomorphic property of the superpotential. That is, there would 
be no way to account for all of the Yukawa terms for both up-type and down-type multiplets 
with a single Higgs doublet. To avoid a massless Higgs state, a mixing term must be 
added to the superpotential. 

From the rules governing the interactions in supersymmetry discussed in the previous 
section, it is easy to see that the terms in ( [82"1) are easily identifiable as fermion masses if the 
Higgses obtain vacuum expectation values (vevs). For example, the first term will contain 
an interaction which we can write as 

1 d 2 W + + 

= -I yeJ fjO(ee c + ete ct ) (84) 
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where it is to be understood that in (53 ) that H\ refers to the scalar component of the Higgs 



Hi and ipL and ip e c represents the fermionic component of the left-handed lepton doublet and 
right-handed singlet respectively. Gauge invariance requires that as defined in fl82|), Hi has 
hypercharge Yh x = — 1 (and Yh 2 = +1). Therefore if the two doublets obtain expectation 
values of the form 

<*> = (?) <*> = (£) (85) 



then (84) contains a term which corresponds to an electron mass term with 

m e = y e vi (86) 

Similar expressions are easily obtained for all of the other massive fermions in the standard 
model. Clearly as there is no v c state in the minimal model, neutrinos remain massless. 
Both Higgs doublets must obtain vacuum values and it is convenient to express their ratio 
parameter of the model, 

tan (3 = — (87) 

Vl 

3.1 The Higgs sector 

Of course if the vevs for Hi and H 2 exist, they must be derivable from the scalar potential 
which in turn is derivable from the superpotential and any soft terms which are included. 
The part of the scalar potential which involves only the Higgs bosons is 

V = \fi\ 2 (H*iHi + H* 2 H 2 ) + lg' 2 (H* 2 H 2 -H*iHi) 2 

o 

+ ^g 2 (A\H*H 2 \ 2 - 2(H*Hi)(H* 2 H 2 ) + {H{Hi) 2 + (H* 2 H 2 ) 2 ) 

+m\HlHi + m\H* 2 H 2 + {B^eijH\H 3 2 + h.c.) (88) 

In ([58]), the first term is a so-called F-term, derived from \(dW / dH{)\ 2 and \(dW / dH 2 )\ 2 
setting all sfermion vevs equal to 0. The next two terms are D-terms, the first a U(l)-D- 
term, recalling that the hypercharges for the Higgses are Yh 1 = — 1 and Yh 2 = 1, and the 
second is an SU(2)-Z)-term, taking T a = <j a /2 where a a are the three Pauli matrices. Finally, 
the last three terms are soft supersymmetry breaking masses mi and m 2 , and the bilinear 
term Bfi. The Higgs doublets can be written as 



w> = UrJ < ff2> = U°> (89) 



and by {H*Hi), we mean tfftf? + Hi*Hi etc. 

The neutral portion of (EST) can be expressed more simply as 



V = ^^(\H i\ 2 -\H 2 \ 2 )\(m 2 i + \fi\ 2 )\H°i\ 2 
8 ^ ' 

+ {m 2 2 + |//| 2 )|# 2 °| 2 + (BfiHiH% + h.c.) (90) 
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For electroweak symmetry breaking, it will be required that either one (or both) of the soft 
masses (m\,m?j>) be negative (as in the standard model). 

In the standard model, the single Higgs doublet leads to one real scalar field, as the other 
three components are eaten by the massive electroweak gauge bosons. In the supersymmetric 
version, the eight components result in 2 real scalars (h,H)\ 1 pseudo-scalar (A); and one 
charged Higgs (i? ); the other three components are eaten. Also as in the standard model, 
one must expand the Higgs doublets about their vevs, and we can express the components 
of the Higgses in terms of the mass eigenstates 

/ v\ + -ks [H cos a — h sin a + iA sin (3] \ 
{Hl)= { ^ H~ sin (3 J 



(H 2 ) 



V2 



H + cos (3 
-j= [H sin a + h cos a + iA cos (3\ 



(91) 



From the two vevs, we can define v 2 = vf + v 2 so that M 2 = \g 2 v 2 as in the standard model. 

In addition, electroweak symmetry breaking places restrictions on the set of parameters 
appearing in the Higgs potential (PO). If we use the two conditions 



dV 







dV 







(92) 



(93) 



(94) 



d\m\ d\m\ 

with a little algebra, we can obtain the following two conditions 

- 2Bn = {ml + m 2 2 + 2fi 2 ) sin 2(3 

and 

2 _ 4 (ml + /i 2 - (ml + fi 2 ) tan 2 (3) 
V = (9 2 + S /2 )(tan 2 - 1) 

From the potential and these two conditions, the masses of the physical scalars can be 
obtained. At the tree level, 

(95) 
(96) 

(97) 
(98) 



m 2 H± 



2 i 2 

m A + m w 



m 



H.h 



m 2 A = ml + m 2 . + 2ji 2 = —Bfx(tan (3 + cot (3) 
1 

" 2 



m A + m z ± J (m\ + m|) 2 — Am 2 A m 2 z cos 2 2(3 



The Higgs mixing angle is defined by 

tan 2a = tan 2(3 



m 2 H + m 2 h 



m , 



TO; 



Notice that these expressions and the above constraints limit the number of free inputs in 
the MSSM. First, from the mass of the pseudoscalar, we see that B[i is not independent and 
can be expressed in terms of m A and tan (3. Furthermore from the conditions (|93| ) and (|9"4]) . 
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we see that if we keep tan/3, we can either either choose tua and [l as free inputs thereby 
determining the two soft masses, rri\ and m 2 , or we can choose the soft masses as inputs, 
and fix rriA and \i by the conditions for electroweak symmetry breaking. Both choices of 
parameter fixing are widely used in the literature. 

The tree level expressions for the Higgs masses make some very definite predictions. 
The charged Higgs is heavier than M\y, and the light Higgs h, is necessarily lighter than 
M z . Note if uncorrected, the MSSM would already be excluded (see discussion on current 
accelerator limits in section 6). However, radiative corrections to the Higgs masses are not 
negligible in the MSSM, particularly for a heavy top mass m t ~ 175 GeV. The leading 
one-loop corrections to m\ depend quartically on m t and can be expressed as [|16|1 



Am 2 



3mf 
Att 2 v 2 



In 



8tt 2 v 2 



2h(ml,ml)+A 2 <7(m|,m|) 



+ 



(99) 



where 2 are the physical masses of the two stop squarks t\$ to be discussed in more detail 
shortly, A t = A t + [l cot (3, (A t is supersymmetry breaking trilinear term associated with the 
top quark Yukawa coupling). The functions h and / are 



a — b \b 



(100) 



Additional corrections to coupling vertices, two-loop corrections and 

renormalization-group resummations have also been computed in the MSSM [P7]. With 
these corrections one can allow 

m h < 130 GeV (101) 

within the MSSM. While certainly higher than the tree level limit of Mz, the limit still 
predicts a relatively light Higgs boson, and allows the MSSM to be experimentally excluded 
(or verified!) at the LHC. 

Recalling the expression for the electron mass in the MSSM flSBD, we can now express the 
Yukawa coupling y e in terms of masses and (tan)/?, 



lie 



V2M W cos 



(102) 



There are similar expressions for the other fermion masses, with the replacement cos/3 
sin (3 for up-type fermions. 



3.2 The sfermions 



We turn next to the question of the sfermion masses [|T^]. As an example, let us consider 
the u mass 2 matrix. Perhaps the most complicated entry in the mass 2 matrix is the L — L 
component. To begin with, there is a soft supersymmetry breaking mass term, m\ 



Q- 



In 
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addition, from the superpotential term, y u H 2 Qu c , we obtain an F-term contribution by 
taking dW/du c = y u H 2 Q. Inserting the vev for H 2) we have in the F-term, 

\dW/du c \ 2 = \y u v 2 u\ 2 = m 2 u \u\ 2 (103) 

This is generally negligible for all but third generation sfermion masses. Next we have the 
D-term contributions. Clearly to generate a u*u term, we need only consider the F-term 
contributions from diagonal generators, i.e., T 3 and Y, that is from 



1 



D 6 = --g^H^ 2 -\H^\ 2 + \u\ 2 + ---\ (104) 

D> = -\g'[\H 2 \ 2 -\H 1 \ 2 + Y Q \u\ 2 + ...} (105) 

where Yq = 2q u — 1 = 1/3 is the quark doublet hypercharge. Once again, inserting vevs for 
Hi and H 2 and keeping only relevant terms, we have for the D-term contribution 

\\D*\ 2 + \\D'\ 2 = \ (g 2 v 2 cos2/3 - g' 2 v 2 cos2(3(2q u - 1)) u*u 

= cos 2 (3 (^-q u sin 2 6 w ^ju*u (106) 

Thus the total contribution to the L — L component of the up-squark mass 2 matrix is 

m 2 Q + m 2 u + M 2 cos 2(3 Q - 5« sin 2 (107) 



Similarly it is easy to see that the R — R component can be found from the above expressions 
by discarding the SU(2)l D-term contribution and recalling that Y u c = 2q u . Then, 

™>l R = m 2 uC +m 2 u + M 2 z cos 2(3 (q u sin 2 W ) (108) 

There are, however, in addition to the above diagonal entries, off-diagonal terms coming 
from a supersymmetry breaking A-term, and an F-term. The A-term is quickly found 
from Aqy u H 2 Qu c when setting the vev for H 2 = v 2 and yields a term Aqm u . The F-term 
contribution comes from dW/dH 2 = [iH.\ + y u Qu c . When inserting the vev and taking 
the square of the latter term, and keeping the relevant mass term, we find for the total 
off-diagonal element 

m l L u R = m u (A Q + ficotp) = m u A Q (109) 

Note that for a down-type sfermion, the definition of A is modified by taking cot j3 — > tan f3. 
Also note that the off-diagonal term is negligible for all but the third generation sfermions. 

Finally to determine the physical sfermion states and their masses we must diagonalize 
the sfermion mass matrix. This mass matrix is easily diagonalized by writing the diagonal 
sfermion eigenstates as 

fi = hcosOf + f R sm6 f , 

f 2 = -f L sm9 f + f R cos9 f . (110) 
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With these conventions we have the diagonalizing angle and mass eigenvalues 



9 



f = sign[-m LR j 



\ — \ arctan 1 2m 2 LR / {m 2 L — m 2 R ) | , m 2 L > m 



R > 



ml<m 2 R , 



| arctan \2m 2 LR J (m| — m 2 R ) \ , 

1,2 - 2 ( m R + m L) T \/( m fi- m i) 2 + 4m l 

Here is chosen so that mi is always lighter that m 2 . Note that in the special case 
mi = m R = m, we have Of = sign[— m| /J? ]('7r/4) and m\ 2 = m 2 =)= 



rn 



l LR 
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3.3 Neutralinos 

There are four new neutral fermions in the MSSM which not only receive mass but mix 
as well. These are the gauge fermion partners of the neutral B and W 3 gauge bosons, 
and the partners of the Higgs. The two gauginos are called the bino, B, and wino, W 3 
respectively. The latter two are the Higgsinos, Hi and if 2- In addition to the supersymmetry 
breaking gaugino mass terms, -\M X BB, and -\M 2 W 3 W 3 , there are supersymmetric mass 
contributions of the type W^ipiipj, giving a mixing term between Hi and H 2 , \[iHiH 2 ^ as 
well as terms of the form g(<j)*T a, ifj)\ a giving the following mass terms after the appropriate 
Higgs vevs have been inserted, -^g'viHiB, —^g'v2H 2 B, —^gviHiW 3 , and -^gv 2 H 2 W 3 . 

These latter terms can be written in a simpler form noting that for example, g'vi/\/2 = 
Mz sin 9y/ cos (3. Thus we can write the neutralino mass matrix as (in the (B,W 3 ,H®,H2) 
basis) 



Mi 




U 1\ 

-M z s 8w cos (3 M z c 8] 
V M z s dw sin (3 -M z c 





M 2 



-M Z Sg w COS (3 
M Z Cg w COS (3 





Mi 



e w cos (3 

e w sin(3 -M z c ew sin f3 -fi 



r z se w sm(3 \ 
M z c dw sin/3 







(112) 



/ 



where sg w = singly and cq w = cos6*iy- The mass eigenstates (a linear combination of the 
four neutralino states) and the mass eigenvalues are found by diagonalizing the mass matrix 
( |112|) . However, by a change of basis involving two new states |19 



S° = Hi sin (3 + H 2 cos (3 

A = - H% cos /3 + H% sin (3 
the mass matrix simplifies and becomes ( in the (B, W 3 , A, S) basis) 

(Mi M z s dw \ 



(113) 
(114) 



Mi 
M 2 

M z s dw -M z c e 




M z s dw 
-M z c ew 

fi sin 2(3 fi cos 2/3 

fi cos 2(3 —fM sin 2(3 ) 



(115) 



In this basis, the eigenstates (which as one can see depend only the the three input mass, 
Mi, M 2 , and \i) can be solved for analytically. 
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Before moving on to discuss the chargino mass matrix, it will be useful for the later 
discussion to identify a few other neutralino states. These are the photino, 

j = W 3 sm6 w + Bcos9 w (116) 

and a symmetric and antisymmetric combination of Higgs bosons, 

H (12) = ^=(H 1 + H 2 ) (117) 
H [ i 2] = : ^(Hi-H 2 ) (118) 



3.4 Charginos 



There are two new charged fermionic states which are the partners of the W gauge bosons 
and the charged Higgs scalars, H ± , which are the charged gauginos, W ± and charged Hig- 
gsinos, H , or collectively charginos. The chargino mass matrix is composed similarly to 
the neutralino mass matrix. The result for the mass term is 



- 5 (w 



H 



f M 2 V^mwsmP^ (~ + ) 



+ h.c. 



(119) 



V \J 7 lm- w cos f3 /i 

Note that unlike the case for neutralinos, two unitary matrices must be constructed to 
diagonalize (|119|) . The result for the mass eigenstates of the two charginos is 

= ~ [Ml + /i 2 + 2M 2 W 

T J (Ml + + 2Mf v y - 4(/iM 2 - ~m\ sin 2f3) 2 1 (120) 



2 2 
Cl ' c 2 



3.5 More Supersymmetry Breaking 

As was noted earlier, supersymmetry breaking can be associated with a positive vacuum 
energy density, V > 0. Clearly from the definition of the scalar potential, this can be 
achieved if either (or both) the F-terms or the D-terms are non-zero. Each of the these two 
possibilities is called F-breaking and D-breaking respectively (for obvious reasons). 



3.5.1 D-Breaking 

One of the simplest mechanisms for the spontaneous breaking of supersymmetry, proposed 
by Fayet and Illiopoulos [f2U| , involves merely adding to the Lagrangian a term proportional 
to D, 

C F i = kD (121) 

It is easy to see by examining ( |73|) that this is possible only for a U(l) gauge symmetry. For 
a U(l), the variation of ( |121| ) under supersymmetry is simply a total derivative. The scalar 
potential is now modified 

V(D) = -i|rf — kD — gig^^D (122) 
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where % is the U(l) charge of the scalar fa. As before, the equations of motion are used to 
eliminate the auxiliary field D to give 

D = -K-g( qi fa*fa) (123) 

So long as the U(l) itself remains unbroken (and the scalar fields fa do not pick up expectation 
values, we can set (fa) = 0, and hence 

(D) = — k (124) 

with 

V=\k 2 (125) 

and supersymmetry is broken. Unfortunately, it is not possible that D-breaking of this type 
occurs on the basis of the known U(l) in the standard model, i.e., U(l)y, as the absence of 
the necessary mass terms in the superpotential would not prevent the known sfermions from 
acquiring vevs. It may be possible that some other U(l) is responsible for supersymmetry 
breaking via .D-terms, but this is most certainly beyond the context of MSSM. 

3.5.2 F-Breaking 

Although F-type breaking also requires going beyond the standard model, it does not do 
so in the gauge sector of the theory. F-type breaking can be achieved relatively easily by 
adding a few gauge singlet chiral multiplets, and one the simplest mechanisms was proposed 



by O'Raifertaigh In one version of this model, we add three chiral supermultiplets, 



A, B, and C, which are coupled through the superpotential 

W = aAB 2 + (3C(B 2 - m 2 ) (126) 

The scalar potential is easily determined from (|126| ) 

V = \F A \ 2 + \F B \ 2 + \F C \ 2 

= \aB 2 \ 2 + \2B(aA + f3C)\ 2 + \(3(B 2 -m 2 )\ 2 (127) 

Clearly, the first and third terms of this potential can not be made to vanish simultaneously, 
and so for example, if B = 0, Fq ^ 0, V > 0, and supersymmetry is broken. 

It is interesting to examine the fermion mass matrix for the above toy model. The mass 
matrix is determined from the superpotential through W^ipiipj and in the (A, B, C) basis 
gives 

/ 2aB \ 

2aB 2(aA + (3C) 2{3B (128) 
V 2(3B / 

The fact that the determinant of this matrix is zero, indicates that there is at least one 
massless fermion state, the Goldstino. 

The existence of the Goldstino as a signal of supersymmetry breaking was already men- 
tioned in the previous section. It is relatively straightforward to see that the Goldstino can 
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be directly constructed from the F- and D-terms which break supersymmetry. Consider the 
mass matrix for a gaugino A a , and chiral fermion ^ 



/ V2gmT a Y\ (uq) 

where we do not assume any supersymmetry breaking gaugino mass. Consider further, the 
fermion 

G = ( (D«)/V2, (Ft) ) (130) 

in the (A, ip) basis. Now from the condition ( |77D and the requirement that we are sitting at 
the minimum of the potential so that 

|K = o <-> g(((j)*}T a y + F 3 W ij = (131) 



we see that the fermion G is massless, that is, it is annihilated by the mass matrix 
The Goldstino state G is physical so long as one or both (D) ^ 0, or (F) ^ 0. This is the 
analog of the Goldstone mechanism for the breaking of global symmetries. 



3.6 R-Parity 

In defining the supersymmetric standard model, and in particular the minimal model or 
MSSM, we have limited the model to contain a minimal field content. That is, the only new 
fields are those which are required by supersymmetry. In effect, this means that other than 
superpartners, only the Higgs sector was enlarged from one doublet to two. However, in 
writing the superpotential flS2"D, we have also made a minimal choice regarding interactions. 
We have limited the types of interactions to include only those required in the standard 
model and its supersymmetric generalization. 

However, even if we stick to the minimal field content, there are several other super- 
potential terms which we can envision adding to fl32|) which are consistent with all of the 
symmetries (in particular the gauge symmetries) of the theory. For example, we could con- 
sider 

W R = \\^ k UL 3 e c k + A"' 7 ' l.,Q ,(11 + Ix'^uffidt + fx H L t H u (132) 
2 2 

In ( |132|) , the terms proportional to A, A', and //, all violate lepton number by one unit. The 
term proportional to A" violates baryon number by one unit. 

Each of the terms in (|132|) predicts new particle interactions and can be to some extent 
constrained by the lack of observed exotic phenomena. However, the combination of terms 
which violate both baryon and lepton number can be disastrous. For example, consider the 
possibility that both A' and A" were non-zero. This would lead to the following diagram 
which mediates proton decay, p — ► e + 7r°, yU + vr°, vn + , vK + etc. Because of the necessary 
antisymmetry of the final two flavor indices in A", there can be no d c exchange in this 
diagram. The rate of proton decay as calculated from this diagram will be enormous due to 
the lack of any suppression by superheavy masses. There is no GUT or Planck scale physics 
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Figure 3: R-parity violating contribution to proton decay. 

which enters in, this is a purely (supersymmetric) standard model interaction. The (inverse) 
rate can be easily estimated to be 

~ 4 
777 

IT 1 ~ lO^eV- 1 (133) 

v m\ 

assuming a supersymmetry breaking scale of m of order 100 GeV. This should be compared 
with current limits to the proton life-time of > 10 63 GeV -1 . 

It is possible to eliminate the unwanted superpotential terms by imposing a discrete 
symmetry on the theory. This symmetry has been called -R-parity |22| , and can be defined 

as 

R = (_l)3B+L+2s ( 134) 

where B, L, and s are the baryon number, lepton number, and spin respectively. With this 
definition, it turns out that all of the known standard model particles have -R-parity +1. 
For example, the electron has B = 0, L = —1, and s = 1/2, the photon as B = L = 
and s = 1. In both cases, R = 1. Similarly, it is clear that all superpartners of the known 
standard model particles have R = —1, since they must have the same value of B and L but 
differ by 1/2 unit of spin. If .R-parity is exactly conserved, then all four superpotential terms 
in ( |132| ) must be absent. But perhaps an even more important consequence of -R-parity is 
the prediction that the lightest supersymmetric particle or LSP is stable. In much the same 
way that baryon number conservation predicts proton stability, -R-parity predicts that the 
lightest R = — 1 state is stable. This makes supersymmetry an extremely interesting theory 
from the astrophysical point of view, as the LSP naturally becomes a viable dark matter 
candidate [^, This will be discussed in detail in the 6th lecture. 

4 The Constrained MSSM and Supergravity 

As a phenomenological model, while the MSSM has all of the ingredients which are necessary, 
plus a large number of testable predictions, it contains far too many parameters to pin down 
a unique theory. Fortunately, there are a great many constraints on these parameters due 
to the possibility of exotic interactions as was the case for additional .R-violating superpo- 
tential terms. The supersymmetry breaking sector of the theory contains a huge number of 
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potentially independent masses. However, complete arbitrariness in the soft sfermion masses 
would be a phenomenological disaster. For example, mixing in the squark sector, would lead 
to a completely unacceptable preponderance of flavor changing neutral currents [24]. 

Fortunately, there are some guiding principles that we can use to relate the various soft 
parameters which not only greatly simplifies the theory, but also removes the phenomeno- 
logical problems as well. Indeed, among the motivations for supersymmetry was a resolution 
of the hierarchy problem ||. We can therefore look to unification (grand unification or even 



unification with gravity) to establish these relations [25]. 

The simplest such hypothesis is to assume that all of the soft supersymmetry breaking 
masses have their origin at some high energy scale, such as the GUT scale. We can further 
assume that these masses obey some unification principle. For example, in the context of 
grand unification, one would expect all of the members of a given GUT multiplet to receive a 
common supersymmetry breaking mass. For example, it would be natural to assume that at 
the GUT scale, all of the gaugino masses were equal, M\ = M 2 = M 3 (the latter is the gluino 
mass) . While one is tempted to make a similar assumption in the case of the sfermion masses 
(and we will do so), it is not as well justified. While one can easily argue that sfermions 
in a given GUT multiplet should obtain a common soft mass, it is not as obvious why all 
R — — 1 scalars should receive a common mass. 

Having made the assumption that the input parameters are fixed at the GUT scale, 
one must still calculate their values at the relevant low energy scale. This is accom- 
plished by "running" the renormalization group equations p6|. Indeed, in standard (non- 
super symmetric) GUTs, the gauge couplings are fixed at the unification scale and run down 
to the electroweak scale. Conversely, one can use the known values of the gauge couplings 
and run them up to determine the unification scale (assuming that the couplings meet at a 
specific renormalization scale). 



4.1 RG evolution 

To check the prospects of unification in detail requires using the two-loop renormalization 
equations 

^ = — ^ + rWk 2 (135) 



where t = \b.(Mq UT /Q 2 ), and the 6, are given by 




x 3 ' 



H 



(136) 
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from gauge bosons, N g matter generations and Nh Higgs doublets, respectively, and at two 
loops 



(137) 



These coefficients depend only on the light particle content of the theory. 

However, using the known inputs at the electroweak scale, one finds @ that the couplings 
of the standard model are not unified at any high energy scale. This is shown in Figure El 
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Figure 4: RG evolution of the inverse gauge couplings in the standard model [0L ITUf. 



In the MSSM, the additional particle content changes the slopes in the RGE evolution 
equations. Including supersymmetric particles, one finds |27| 
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In this case, it is either a coincidence, or it is rather remarkable that the RG evolution 
is altered in just such a way as to make the MSSM consistent with unification. The MSSM 
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Figure 5: RG evolution of the inverse gauge couplings in the MSSM J|, [10| . 



evolution is shown in Figure [5] below. For many, this concordance of the gauge couplings 
and GUTs offers strong motivation for considering supersymmetry. 

As was noted earlier, most of the parameters in the MSSM are also subject to RG 
evolution. For example, all of the Yukawa couplings also run. Here, I just list the RG 
equation for the top quark Yukawa, 

da t a t / 16 13 \ 

— — = — — «3 + 3«2 H oil — oa t — ah + • ■ ■ ] (140) 

dt 4vr V 3 15 J K J 

where a t = yf/Air. This is the leading part of the 1-loop correction. For a more complete 
list of these equations see ||28|| . These expressions are also known to higher order ||29|| . 



Note that the scaling of the supersymmetric couplings are all proportional to the couplings 
themselves. That means that if the coupling is not present at the tree level, it will not be 
generated by radiative corrections either. This is a general consequence of supersymmetric 



nonrenormalization theorems 30 



The supersymmetry breaking mass parameters also run. Starting with the gaugino 
masses, we have 

^ = -ka.Mi/Arc (141) 

Assuming a common gaugino mass, at the GUT scale as was discussed earlier, these 
equations are easily solved in terms of the fine structure constants, 

ai{M GUT ) 

This implies that 

Mi M-2 M 3 , 
9i 92 9i 
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(Actually, in a GUT, one must modify the relation due to the difference between the U(l) 
factors in the GUT and the standard model, so that we have Mi = |— Mo.) 

Finally, we have the running of the remaining mass parameters. A few examples are: 

df^ 2 3u 2 / 1 \ 

a 2 + -at - a t - a b + ■ ■ ■ (144) 



dt A-k \ 5 

= l ^Ml + ... (145) 

dt 5 47r 1 v ' 

dAi 1/16 13 \ , , 

~~df = 47 (l° 3M3 + 3c * 2M2 + Y^" lMl _ 6a * A * + ' " ) ( 146 ) 

^ = £■ (a 2 M 2 + i ai Mi - 3a*A + • • •) (147) 
4.2 The Constrained MSSM 

As the name implies, the constrained MSSM or CMSSM, is a subset of the possible parameter 



sets in the MSSM. In the CMSSM [[H], we try to make as many reasonable and well 
motivated assumptions as possible. To begin with gaugino mass unification is assumed. 
(This is actually a common assumption in the MSSM as well). Furthermore soft scalar 
mass unification or universality is also assumed. This implies that all soft scalar masses are 
assumed equal at the GUT input scale, so that 

m 2 (M GUT )=m 2 (148) 

This condition is applied not only to the sfermion masses, but also to the soft Higgs masses, 
ml 2 as we h. By virtue of the conditions (0) and (|94]), we see that in the CMSSM, /i, and 
Bji, (or m\), are no longer free parameters since these conditions amount to m^(m 2 , m§, fi 2 ) 
and v 2 ((m1, m?,, /i 2 ). Thus we are either free to pick rriA, A* as f ree parameters (this fixes m^, 
though we are usually not interested in those quantities) as in the MSSM, or choose m 12 
(say at the GUT scale) and and /i become predictions of the model. Universality of the 
soft trilinears, A iy is also assumed. 

In the CMSSM therefore, we have only the following free input parameters: mi/2, mo, 
tan/5, A , and the sign of fi. We could of course choose phases for some these parameters. 
In the MSSM and CMSSM, there are two physical phases which can be non-vanishing, 9^, 
and 9 a- If non-zero, they lead to various CP violating effects such as inducing electric 
dipole moments in the neutron and electron. For some references regarding these phases see 



|33| , |34| , p5| , but we will not discuss them further in these lectures. 

In the figure below, an example of the running of the mass parameters in the CMSSM 
is shown. Here, we have chosen = 250 GeV, m = 100 GeV, tan/? = 3, A = 0, and 
fi < 0. Indeed, it is rather amazing that from so few input parameters, all of the masses 
of the supersymmetric particles can be determined. The characteristic features that one 
sees in the figure, are for example, that the colored sparticles are typically the heaviest in 
the spectrum. This is due to the large positive correction to the masses due to a% in the 
RGE's. Also, one finds that the B, is typically the lightest sparticle. But most importantly, 
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Figure 6: RG evolution of the mass parameters in the CMSSM. I thank Toby Falk for 
providing this figure. 

notice that one of the Higgs mass 2 , goes negative triggering electroweak symmetry breaking 
| 3lfl . (The negative sign in the figure refers to the sign of the mass 2 , even though it is the 



mass of the sparticles which are depicted.) In the Table below, I list some of the resultant 
electroweak scale masses, for the choice of parameters used in the figure. 

4.3 Supergravity 

Up until now, we have only considered global supersymmetry. Recall, our generator of 
supersymmetry transformations, the spinor £. In all of the derivations of the transformation 
and invariance properties in the first two sections, we had implicitly assumed that <9 M £ = 0. 
By allowing £(a;) and <9 M £(x) ^ 0, we obtain local supersymmetry or supergravity [p6| . It is 



well beyond the means of these lectures to give a thorough treatment of local supersymmetry. 
We will therefore have to content ourself with some general remarks from which we can 
glimpse at some features for which we can expect will have some phenomenological relevance. 

First, it is important to recognize that our original Lagrangian for the Wess-Zumino 
model involving a single noninteracting chiral multiplet will no longer be invariant under local 
supersymmetry transformations. New terms, proportional to d^^x) must be canceled. In 
analogy with gauge, theories which contain similar terms and are canceled by the introduction 
of vector bosons, here the terms must be canceled by introducing a new spin 3/2 particle 
called the gravitino. The supersymmetry transformation property of the gravitino must be 

« oc d,C (149) 
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Table 1: Physical mass parameters for 7711/2 = 250 GeV, mo = 100 GeV, tan/3 = 3, Aq = 0, 
and /i < 0. (These are related but not equal to the running mass parameters shown in Figure 

B) 



particle 


mass 


parameter 


value 


m i 


203 


n 


-391 


m e -c 


144 


M 1 


100 




190 


M 2 


193 


rriq 


412-554 


M 3 


726 


m xi 


104 


a 3 (M z ) 


.123 


m xf 
m h 


203 


A's 


163-878 


93 






m A 


466 







Notice that the gravitino carries both a gauge and a spinor index. The gravitino is part of 
an N = 1 multiplet which contains the spin two graviton. In unbroken supergravity, the 
massless gravitino has two spin components (± 3/2) to match the two polarization states of 
the graviton. 

The scalar potential is now determined by a analytic function of the scalar fields, called 
the Kahler potential, K((f),<f)*). The Kahler potential can be thought of as a metric in field 
space, 

Ckin = Kfi^d^* (150) 

where K l = dK/dipi and Ki = dK/dcff* . In what is known as minimal N — 1 supergravity, 
the Kahler potential is given by 

K = k 2 ^* +\ii(k 6 \W\ 2 ) (151) 
where W((j)) is the superpotential, kT 1 = Mp/y8n and the Planck mass is Mp = 1.2 x 10 19 



GeV. The scalar potential (neglecting any gauge contributions) is J37 



v{<t>A* 



K -4 

en 



(152) 



For minimal supergravity, we have K % = K 2 <p l +W l /W } Ki = K 2 (pi + W* /W* , and (K 
6j/K 2 . Thus the resulting scalar potential is 



-IV 



\W* + 6 % W\ 2 - 3k?\W\ 



(153) 



As we will now see, one of the primary motivations for the CMSSM, and scalar mass 
universality comes from the simplest model for local supersymmetry breaking. The model 
|38| involves one additional chiral multiplet z, (above the normal matter fields (pi). Let 
us consider therefore, a superpotential which is separable in the so-called Polonyi field and 
matter so that 

W(z,</> i ) = f{z)+g(<l> i ) (154) 
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and in particular let us choose 

f(z)=fx(z + f3) (155) 

and for reasons to be clear shortly, (3 = 2 — v3. I will from now on work in units such that 
k = 1. If we ignore for the moment the matter fields 0, the potential for z becomes 



V{z,z*) 



e zz *fi 2 



\l + z*(z + (3)\ 2 -3\(z + p)\' 



(156) 



It is not difficult to verify that with the above choice of /3, the minimum of V occurs at 
(z) = V3- 1, with V((z)) = 0. 

Note that by expanding this potential, one can define two real scalar fields A and B, with 
mass eigenvalues, 



2V3 



where the gravitino mass is 



m 3/2 

m 3 / 2 ~- 



m 



B 



2(2 - V3) 



rn. 



3/2 



e K/2 
2 1 ™2 



e 2 ^ 



(157) 



(158) 



Note also that there is a mass relation, m A + m B = 4m^ 2 , which is a guaranteed consequence 
of supertrace formulae in supergravity ||37|| . Had we considered the fermionic sector of the 
theory, we would now find that the massive gravitino has four helicity states ±1/2 and ±3/2. 
The "longitudinal" states arise from the disappearance of the goldstino (the fermionic partner 
of z in this case) by the superHiggs mechanism, again in analogy with the spontaneous 
breakdown of a gauge symmetry |37], . 

We next consider the matter potential from eqs. ( |154|) and (|153|) . In full, this takes the 



form 40 



V 



,(kl 2 +l</>l 2 ) 



dg 



!— 

dz 



+ z*(f(z)+g(cf>))\' 



(159) 



Here again, I have left out the explicit powers of Mp. Expanding this expression, and at the 
same time dropping terms which are suppressed by inverse powers of the Planck scale (this 
can be done by simply dropping terms of mass dimension greater than four), we have, after 
inserting the vev for z ||40|| , 



V 



|//±(V3-l)Gu±s(0))f 



+ 



dg 



±0*(/i±^(0))| 2 -3|/i±^(0)p 



,(4-2^3) 



■V3a*W) + ^W) + ||tI 



,(4-2 



+m 3/2 e 



(2-^3), 



dg 



— V3g + h.c.)) + m 2 



3/2S 



(160) 
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This last expression deserves some discussion. First, up to an overall rescaling of the 
superpotential, g — > e^~~ 2 g, the first term is the ordinary F-term part of the scalar poten- 
tial of global supersymmetry. The next term, proportional to 772,3/2 represents a universal 
trilinear A-term. This can be seen by noting that J24>dg/d<p = 3g, so that in this model 
of supersymmetry breaking, A = (3 — V3)m 3 / 2 . Note that if the superpotential contains 
bilinear terms, we would find B = (2 — \f2t)m^/2- The last term represents a universal scalar 
mass of the type advocated in the CMSSM, with tUq = m 2 , 2 . The generation of such soft 



terms is a rather generic property of low energy supergravity models [41 



Before concluding this section, it worth noting one other class of supergravity models, 



namely the so-called no-scale supergravity model J|2] . No-scale supergravity is based on the 
Kahler potential of the form 

K — — \n{S + S*) - 3 ln(T + T* — faf*) + In \W\ 2 (161) 

where the S and T fields are related to the dilaton and moduli fields in string theory [f43 |. 



If only the T field is kept in ( |161[ ), the resulting scalar potential is exactly flat, i.e., V — 
identically. In such a model, the gravitino mass is undetermined at the tree level, and up 
to some field redefinitions, there is a surviving global supersymmetry. No-scale supergravity 
has been used heavily in constructing supergravity models in which all mass scales below 
the Planck scale are determined radiatively pi|,[4"5 . 



5 Cosmology 

Supersymmetry has had a broad impact on cosmology. In these last two lectures, I will try to 
highlight these. In this lecture, I will briefly review the problems induced by supersymmetry, 
such as the Polonyi or moduli problem, and the gravitino problem. I will also discuss the 
question of cosmological inflation in the context of supersymmetry. Finally, I will describe 
a mechanism for baryogenesis which is purely supersymmetric in origin. I will leave the 
question of dark matter and the accelerator constraints to the last lecture. 

Before proceeding to the problems, it will be useful to establish some of the basic quan- 
tities and equations in cosmology. The standard big bang model assumes homogeneity and 
isotropy, so that space-time can be described by the Friedmann-Robertson- Walker metric 
which in co-moving coordinates is given by 



d s 2 = -dt 2 + R 2 (t) 



- dT ] 2 , + r 2 (^ 2 + sinW) 
1 — kr 2 ) v ' 



(162) 



where R(t) is the cosmological scale factor and k is the three-space curvature constant 
(k = 0, +1,— 1 for a spatially flat, closed or open Universe), k and R are the only two 
quantities in the metric which distinguish it from flat Minkowski space. It is also common 
to assume the perfect fluid form for the energy-momentum tensor 

Tp, = V9iiv + (P + p)un u " ( 163 ) 
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where is the space-time metric described by ( |162| ), p is the isotropic pressure, p is the 
energy density and w M = (1, 0, 0, 0) is the velocity vector for the isotropic fluid. Einstein's 
equation yield the Friedmann equation, 

H 2= (R\ = I K 2 1, I A (164) 
\R 3 H R 2 3 K J 



and 



|) =^A-^ 2 (p + 3p) (165) 



where A is the cosmological constant, or equivalently from T^ v '. v = 

p = -3H(p + p) (166) 

These equations form the basis of the standard big bang model. 

If our Lagrangian contains scalar fields, then from the scalar field contribution to the 
energy-momentum tensor 

Tfa, = d^d u <P - l -g^d p <pd p <P - g^V(<P) (167) 
we can identify the energy density and pressure due to a scalar 0, 

p= l -^+ l -R-\t)iy<Pf + V(<P) (168) 

p = ^ 2 -i J R- 2 (t)(V0) 2 -y(0) (169) 
I 6 

In addition, we have the equation of motion, 

dV 

+ 3#0 + — = (170) 

ocp 

Finally, I remind the reader that in the early radiation dominated Universe, the energy 
density (as well as the Hubble parameter) is determined by the temperature, 



p = N-T i H=J—kT 2 (171) 
H 30 V 90 K J 

The critical energy density (corresponding to k — 1, is 

p c = 3H\- 2 = 1.88 x 10 - 29 gcirT 3 /io (172) 

where the scaled Hubble parameter is h = i/ /100km Mpc _1 s _1 . The cosmological density 
parameter is defined as Q = p/p c . 



38 



5.1 The Polony i Problem 

The Polonyi problem, although based on the specific choice for breaking supergravity dis- 
cussed in the previous lecture (eq. |155|) , is a generic problem in broken supergravity models. 



In fact, the problem is compounded in string theories, where there are in general many such 
fields called moduli. Here, attention will be restricted to the simple example of the Polonyi 
potential. 

The potential in eq. ( |156| ) has the property that at its minimum occurs at (z) = — 
1)M, where M = k,^ 1 is the reduced Planck mass. Recall that the constant (3 was chosen 
so that at the minimum, V((z)) = 0. In contrast to the value of the expectation value, the 
curvature of the potential at the minimum, is m 2 ~ p 2 , which as argued earlier is related 
to the gravitino mass and therefore must be of order the weak scale. In addition the value 
of the potential at the origin is of order V(0) ~ p 2 M 2 , i.e., an intermediate scale. Thus, we 
have a long and very flat potential. 

Without too much difficulty, it is straightforward to show that such a potential can be 
very problematic cosmologically j46|. The evolution of the Polonyi field z, is governed by 
eq. (|170|) with potential ( |156| ). There is no reason to expect that the field z is initially at its 
minimum. This is particularly true if there was a prior inflationary epoch, since quantum 
fluctuations for such a light field would be large, displacing the field by an amount of order 
M from its minimum. If z ^ (z), the evolution of z must be traced. When the Hubble 
parameter H > p , z is approximately constant. That is, the potential term (proportional to 
p 2 ) can be neglected. At later times, as H drops, z begins to oscillate about the minimum 
when H < p. Generally, oscillations begin when H ~ m z ~ /i as can be seen from the 
solution for the evolution of a non-interacting massive field with V = m 2 z z 2 /2. This solution 
would take the form of z ~ sm(m z t)/t with H = 2/3t. 

At the time that the ^-oscillations begin, the Universe becomes dominated by the po- 
tential V(z), since H 2 ~ p/M 2 . Therefore all other contributions to p will redshift away, 
leaving the potential as the dominant component to the energy density. Since the oscillations 
evolve as non-relativistic matter (recall that in the above solution for z, H = 2/3t as in a 
matter dominated Universe). As the Universe evolves, we can express the energy density 
as p ~ p 2 M 2 (R z /R) 3 , where R z is the value of the scale factor when the oscillations begin. 
Oscillations continue, until the ^-fields can decay. Since they are only gravitationally coupled 
to matter, their decay rate is given by T z ~ p 3 /M 2 . Therefore oscillations continue until 
H ~ Y z or when R = Rd z ~ (M//i) 4//3 . The energy density at this time is only p 6 /M 2 . Even 
if the the thermalization of the decay products occurs rapidly, the Universe reheats only to 
a temperature of order T R ~ p 1/4 ~ p 3/2 /M 1/2 . For p ~ 100 GeV, we have T R ~ 100 keV! 
There are two grave problems with this scenario. The first is is that big bang nucleosynthesis 
would have taken place during the oscillations which is characteristic of a matter dominated 
expansion rather than a radiation dominated one. Because of the difference in the expansion 
rate the abundances of the light elements would be greatly altered (see e.g. |]47||). Even more 
problematic is the entropy release due to the decay of these oscillations. The entropy increase 
f46| is related to the ratio of the reheat temperature to the temperature of the radiation in 



the Universe when the oscillations decay, Td ~ Ti(R z / Rd z ) where Tj is the temperature when 
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oscillations began Tj ~ (/iM) 1//2 . Therefore, the entropy increase is given by 

Sf/Si ~ (T^/Trf) 3 ~ (M//i) ~ 10 16 (173) 

This is far too much to understand the present value of the baryon-to-entropy ratio, of order 
10~ n — 1CT 10 as required by nucleosynthesis and allowed by baryosynthesis. That is, even 
if a baryon asymmetry of order one could be produced, a dilution by a factor of 10 16 could 
not be accommodated. 



5.2 The Gravitino Problem 

Another problem which results from the breaking of supergravity is the gravitino problem 
f48fl . If gravitinos are present with equilibrium number densities, we can write their energy 
density as 



P3/2 = ^3/2^3/2 = m 3/2 ^-^T~) T 3/2 ( 174 ) 

where today one expects that the gravitino temperature T 3 / 2 is reduced relative to the 
photon temperature due to the annihilations of particles dating back to the Planck time 
P9| . Typically one can expect Y 3 / 2 = (T 3 / 2 /T 7 ) 3 ~ 10~ 2 . Then for f2 3 / 2 /i 2 < 1, we obtain 
the limit that m 3 / 2 ^ 1 keV. 

Of course, the above mass limits assumes a stable gravitino, the problem persists however, 
even if the gravitino decays, since its gravitational decay rate is very slow. Gravitinos decay 
when their decay rate, T 3 / 2 ~ mL/Mp, becomes comparable to the expansion rate of the 

Universe (which becomes dominated by the mass density of gravitinos), H ~ rn^^T^ /Mp. 

Thus decays occur at T d ~ m^l/Mp. After the decay, the Universe is "reheated" to a 
temperature 

T R * p(T d )V* ~ m%/M P 12 (175) 

As in the case of the decay of the Polonyi fields, the Universe must reheat sufficiently so that 
big bang nucleosynthesis occurs in a standard radiation dominated Universe. For Tr > 1 
MeV, we must require m 3 / 2 > 20 TeV. This large value threatens the solution of the hierarchy 
problem. In addition, one must still be concerned about the dilution of the baryon-to-entropy 
ratio f50fl , in this case by a factor A = {T R /T D f ~ Y(M P /m 3/2 )^ 2 . Dilution may not be a 



problem if the baryon-to-entropy ratio is initially large. 

Inflation (discussed below) could alleviate the gravitino problem by diluting the gravitino 
abundance to safe levels [|50| . If gravitinos satisfy the noninflationary bounds, then their 
reproduction after inflation is never a problem. For gravitinos with mass of order 100 GeV, 
dilution without over-regeneration will also solve the problem, but there are several factors 
one must contend with in order to be cosmologically safe. Gravitino decay products can 
also upset the successful predictions of Big Bang nucleosynthesis, and decays into LSPs (if 
R-parity is conserved) can also yield too large a mass density in the now-decoupled LSPs 
T9| . For unstable gravitinos, the most restrictive bound on their number density comes form 



the photo-destruction of the light elements produced during nucleosynthesis |>T| 

n 3 / 2 /n 7 < 10- 13 (100GeV/m 3/2 ) (176) 
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for lifetimes > 10 4 sec. Gravitinos are regenerated after inflation and one can estimate 



ia 50, 51 



n V2 /n, ~ (r/#)(T 3/2 /T 7 ) 3 ~ aN(T R )(T R /M P )(T 3/2 /T 7 ) 3 (177) 

where T ~ aN(T R )(Tp/Mp) is the production rate of gravitinos. Combining these last two 
equations one can derive bounds on T R 

T R <Ax 10 9 GeV(100 GeV/m 3/2 ) (178) 



using a more recent calculation of the gravitino regeneration rate p2 \. A slightly stronger 



bound (by an order of magnitude in T R ) was found in |53] 



5.3 Inflation 

It would be impossible in the context of these lectures to give any kind of comprehensive 
review of inflation whether supersymmetric or not. I refer the reader to several reviews 
|54]| . Here I will mention only the most salient features of inflation as it connects with 



sup er symmet ry. 

Supersymmetry was first introduced [55] in inflationary models as a means to cure some 
of the problems associated with the fine-tuning of new inflation New inflationary 

models based on a Coleman- Weinberg type of SU(5) breaking produced density fluctuations 



[57J with magnitude <5p/p~O(10 2 ) rather than <5p/p ~ 10 5 as needed to remain consistent 



with microwave background anisotropics. Other more technical problems [pq] concerning slow 
rollover and the effects of quantum fluctuations also passed doom on this original model. 

The problems associated with new inflation, had to with the interactions of the scalar field 
driving inflation, namely the SU(5) adjoint. One cure is to (nearly) completely decouple the 
field driving inflation, the inflaton, from the gauge sector. As gravity becomes the primary 
interaction to be associated with the inflaton it seemed only natural to take all scales to 
be the Planck scale [0. Supersymmetry was also employed to give flatter potentials and 



acceptable density perturbations [p9| . These models were then placed in the context of N=l 
supergravity p0|, |5"1"| ]. 

The simplest such model for the inflaton 77, is based on a superpotential of the form 

W{rj) = p 2 (l - r]/Mp) 2 M P (179) 



or 



W{ri) = fi 2 (v ~ V 4 /4:M P 3 ) (180) 



where Eq.( |179| )pT| is to be used in minimal supergravity while Eq. ( |180| )p^[ is to be used in 



no-scale supergravity. Of course the real goal is to determine the identity of the inflaton. 
Presumably string theory holds the answer to this question, but a fully string theoretic 



inflationary model has yet to be realized [63 



For the remainder of the discussion, it will be sufficient to consider only a generic model 
of inflation whose potential is of the form: 

V(rj) = n A P{rj) (181) 
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where 77 is the scalar field driving inflation, the inflaton, \i is an as yet unspecified mass 
parameter, and P{rj) is a function of rj which possesses the features necessary for inflation, 
but contains no small parameters, i.e., where all of the couplings in P are 0(1) but may 
contain non-renormalizable terms. 

The requirements for successful inflation can be expressed in terms of two conditions: 
1) enough inflation; 

d 2 V 3H 2 8vrU(0) 

\r 1 ~r ll ±H< — — — 2 



drf 
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2) density perturbations of the right magnitude [57| 
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H 2 



10tt 3 / 2 ?7 



0(100) 
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given here for scales which "re-enter" the horizon during the matter dominated era. For 
large scale fluctuations of the type measured by COBE[54], we can use Eq. ( |183| ) to fix the 
inflationary scale \i p5 |: 



M 2 



few x 10" 



(184) 

For example, 
10 ■.]!,,. The 

duration of inflation is r ~ Mp//i 4 , and the number of e-foldings of expansion is Hr ~ 
87r(M|,//i 2 ) ~ 10 9 . If the inflaton decay rate goes as T ~ m^/Mp ~ fi 6 /Mp, the universe 

' n M P ~ 10 8 GeU. However, it was 



Fixing (n 2 / Mp) has immediate general consequences for inflationpB|. 
the Hubble parameter during inflation, H 2 ~ (87r/3)(/z 4 /Mp) so that H ~ 



recovers at a temperature T R ~ (TM P 



iV2 



10" 



noted in |6^] that in fact the Universe is not immediately thermalized subsequent to inflaton 
decays, and the process of thermalization actually leads to a smaller reheating temperature, 



T R ~ a 2 ^/M 2 P ~ l^GeV 



(185) 



where a 2 ~ 10 -3 characterizes the strength of the interactions leading to thermalization. 
This low reheating temperature is certainly safe with regards to the gravitino limit (|178|) 
discussed above. 



5.4 Baryogenesis 

The production of a net baryon asymmetry requires baryon number violating interactions, 
C and CP violation and a departure from thermal equilibrium^]. The first two of these 
ingredients are contained in GUTs, the third can be realized in an expanding universe where 
it is not uncommon that interactions come in and out of equilibrium. 



In the original and simplest model of baryogenesis ||68|| , a GUT gauge or Higgs boson 
decays out of equilibrium producing a net baryon asymmetry. While the masses of the 
gauge bosons is fixed to the rather high GUT scale 10 15 ~ 16 GeV, the masses of the triplets 
could be lighter O(10 10 ) GeV and still remain compatible with proton decay because of 
the Yukawa suppression in the proton decay rate when mediated by a heavy Higgs. This 
reduced mass allows the simple out-out-equilibrium decay scenario to proceed after inflation 
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so long as the Higgs is among the inflaton decay products p9|. From the arguments above, 
an inflaton mass of 10 11 GeV is sufficient to realize this mechanism. Higgs decays in this 
mechanism would be well out of equilibrium as at reheating T <C ran and nn ~ n 
case, the baryon asymmetry is given simply by 



7- 



In this 
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M, 



10" 4 e 



(186) 



where e is the CP violation in the decay, Tr is the reheat temperature after inflation, and I 



have substituted for n v = p ri /m ri ~ T Mp /m v . 



In a supersymmetric grand unified SU(5) theory, the superpotential Fy must be expressed 
in terms of SU(5) multiplets 



Y 



h d H 2 5 10 + huH.! 10 10 
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where 10, 5, H\ and H 2 are chiral supermultiplets for the 10 and 5-plets of SU(5) matter fields 
and the Higgs 5 and 5 multiplets respectively There are now new dimension 5 operators 
25|, [?T| which violate baryon number and lead to proton decay as shown in Figure |^. The 



first of these diagrams leads to effective dimension 5 Lagrangian terms such as 



'(5) 
-eff 



-w {qqql) 



H 



and the resulting dimension 6 operator for proton decay 
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(188) 



(189) 



As a result of diagrams such as these, the proton decay rate scales as V ~ h 4 g A /Mf I M^ ; 
where Mh is the triplet mass, and M G is a typical gaugino mass of order < 1 TeV. This rate 
however is much too large unless M H > 10 16 GeV. 
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Figure 7: Dimension 5 and induced dimension 6 graphs violating baryon number. 

It is however possible to have a lighter (O(10 10 — 10 11 ) GeV) Higgs triplet needed for 
baryogenesis in the out-of-equilibrium decay scenario with inflation. One needs two pairs of 
Higgs five-plets (Hi, H 2 and H[, H' 2 which is anyway necessary to have sufficient C and CP 
violation in the decays. By coupling one pair (H 2 and H[) only to the third generation of 
fermions via [[f^| 

aHilOlO + feHilOalOs + cHalOaSa + c/H'alOB (190) 
proton decay can not be induced by the dimension five operators. 
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5.4.1 The Affleck-Dine Mechanism 



Another mechanism for generating the cosmological baryon asymmetry is the decay of scalar 
condensates as first proposed by Affleck and Dine[fn|. This mechanism is truly a product of 
super symmetry. It is straightforward though tedious to show that there are many directions 
in field space such that the scalar potential given in eq. ([79]) vanishes identically when SUSY 
is unbroken. That is, with a particular assignment of scalar vacuum expectation values, 
V = in both the F— and D — terms. An example of such a direction is 



?i<t> 



\Jv 2 + a 2 



(191) 



where a, v are arbitrary complex vacuum expectation values. SUSY breaking lifts this de- 
generacy so that 

V ~ m 2 2 (192) 



where m is the SUSY breaking scale and <fi is the direction in field space corresponding to 
the flat direction. For large initial values of <f), <p 



4> ~ M guU a large baryon asymmetry can 
be generated [fT3|, 7^| . This requires the presence of baryon number violating operators such 
as O = qqql such that (O) ^ 0. The decay of these condensates through such an operator 
can lead to a net baryon asymmetry. 

In a supersymmetric gut, as we have seen above, there are precisely these types of op- 
erators. In Figure [8], a 4-scalar diagram involving the fields of the flat direction ( |191[ ) is 
shown. Again, G is a (light) gaugino, and X is a superheavy gaugino. The two supersym- 
metry breaking insertions are of order m, so that the diagram produces an effective quartic 
coupling of order m 2 /{cf) 2 + Mj^). 
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Figure 8: Baryon number violating diagram involving flat direction fields. 



The baryon asymmetry is computed by tracking the evolution of the sfermion condensate, 
which is determined by 

<j> + 3H<j> = -m 2 (j) (193) 
To see how this works, it is instructive to consider a toy model with potential |74 



U(0,0*)=m 2 #* + 2^^ 4 -^ 4 ] 



(194) 
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The equation of motion becomes 



01 + 3H4>t = -m 2 0! + 3A0 2 2 _ (195) 

2 + 3#0 2 = -m 2 2 - 3A0 2 0i + A0? (196) 

with = (0i + z0 2 ) / \/2. Initially, when the expansion rate of the Universe, H, is large, we 
can neglect and fh. As one can see from (|194j) the flat direction lies along ~ 0i ~ O 



with 02 ~ 0. In this case, 0i ~ and 02 — Since the baryon density can be written as 

n B = jo = ^(0102 — 0201 ) — 67?0o; by generating some motion in the imaginary direction, 
we have generated a net baryon density. 

When H has fallen to order fh (when ~ m), 0i begins to oscillate about the origin 
with 0i ~ O sm(fht) / fht . At this point the baryon number generated is conserved and the 
baryon density, n B falls as R~ 3 . Thus, 

n B ~ — 2 2 oc R~ 3 (197) 
m 

and relative to the number density of 0's (n^ = p^/fh = fhcf) 2 ) 

n B A0 2 



fh 2 



(198) 



If it is assumed that the energy density of the Universe is dominated by 0, then the 
oscillations will cease, when 

om 

or when the amplitude of oscillations has dropped to 0^ ~ (Mpfh 2 ) 1 ^ 3 . Note that the decay 
rate is suppressed as fields coupled directly to gain masses oc 0. It is now straightforward 
to compute the baryon to entropy ratio, 

-^^%.M(^y /6 (200) 

s p 3 J* fh 5 / 2 ^ 2 m 2 \m J 
and after inserting the quartic coupling 

*»- ^ ' M ^'° (201) 



s {M\ + 02) V fh 
which could be 0(1). 

In the context of inflation, a couple of significant changes to the scenario take place. 
First, it is more likely that the energy density is dominated by the inflaton rather than 
the sfermion condensate. The sequence of events leading to a baryon asymmetry is then 



as follows pS |: After inflation, oscillations of the inflaton begin at R = R n when H 



and oscillations of the sfermions begin at R = R<f, when H ~ fh. If the Universe 
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is inflaton dominated 
can relate R n and R, 

r 



p l J 2 and 



'r] """" Pv 

As discussed earlier 



H ~ m^R^/Rf/ 2 since # 
-nu R^ ~ (m v /rh) 2/3 R v . 
,, - m^/M| = H or when R = R dri ~ (M p /m v ) 4 / 3 R v . 

dominated by the relativistic decay products of the inflaton, p rr) = m 2 ^Mp ^ (R^/RY and 

H = m l^M^iR^/R) 2 . Sfermion decays still occur when = H which now corresponds 

to a value of the scale factor R^ = (ml/ 11 '(j) 2 /* 'Mp 15 ' / 'mjR^. The final baryon asymmetry in 
the Affleck-Dine scenario with inflation becomes 



r] ~ R Thus one 
inflatons decay when 
The Universe then becomes 

. _ ™2/ 
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(202) 



for m ~ (10~ 17 - 10- 16 )M P , and M x ~ (10~ 4 - 10- 3 )Mp and m v ~ (10~ 8 - 10- ? )M P . 

When combined with inflation, it is important to verify that the AD flat directions remain 
flat. In general, during inflation, supersymmetry is broken. The gravitino mass is related to 
the vacuum energy and m 2 , 2 ~ V/Mp ~ H 2 , thus lifting the flat directions and potentially 
preventing the realization of the AD scenario as argued in | 75[| . To see this, recall the minimal 
supergravity model defined in eqs. (|151|) - (|154|) . Recall also, the last term in eq. ( |160| ), 
which gives a mass to all scalars (in the matter sector), including flat directions of order 
the gravitino mass which during inflation is large. This argument can be generalized to 
supergravity models with non-minimal Kahler potentials. 

However, in no-scale supergravity models, or more generally in models which possess a 
Heisenberg symmetry |76[1 , the Kahler potential can be written as (cf. eq. ( |161| )) 



G = f{z + z* 



Now, one can write 
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(204) 



It is important to notice that the cross term |<^*W| 2 has disappeared in the scalar potential. 
Because of the absence of the cross term, flat directions remain flat even during inflation 
77 1 . The no-scale model corresponds to / = —3 In 77, f' 2 = 3f" and the first term in (|204|) 



vanishes. The potential then takes the form 



V 



e*'\WA 



(205) 



which is positive definite. The requirement that the vacuum energy vanishes implies (Wj) = 
(g a ) = at the minimum. As a consequence 7] is undetermined and so is the gravitino mass 

^3/2 (?7)- 



The above argument is only valid at the tree level. An explicit one-loop calculation |78 
shows that the effective potential along the flat direction has the form 
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where A is the cutoff of the effective supergravity theory, and has a minimum around ~ 
0.5A. Thus, 0o ~ M P will be generated and in this case the subsequent sfermion oscillations 
will dominate the energy density and a baryon asymmetry will result which is independent 



of inflationary parameters as originally discussed in |73], |74J and will produce %/s ~ 0(1) 



Thus we are left with the problem that the baryon asymmetry in no-scale type models is too 



large [77 



In ||80|1 , several possible solutions were presented to dilute the baryon asymmetry. These 
included 1) entropy production from moduli decay, 2) the presence of non-renormalizable 
interactions, and 3) electroweak effects. Moduli decay in this context, turns out to be 
insufficient to bring an initial asymmetry of order Ub/s ~ 1 down to acceptable levels. 
However, as a by-product one can show that there is no moduli problem either. In contrast, 
adding non-renormalizable Planck scale operators of the form <fi 2n ~ 2 / Mp n ~ G leads to a smaller 
initial value for <p and hence a smaller value for ub/s. For dimension 6 operators (n = 4), a 
baryon asymmetry of order n^/s ~ 1CT 10 is produced. Finally, another possible suppression 
mechanism is to employ the smallness of the fermion masses. The baryon asymmetry is 
known to be wiped out if the net B — L asymmetry vanishes because of the sphaleron 
transitions at high temperature. However, Kuzmin, Rubakov and Shaposhnikov |81[] pointed 
out that this erasure can be partially circumvented if the individual (B — 3Lj) asymmetries, 
where i = 1, 2, 3 refers to three generations, do not vanish even when the total asymmetry 
vanishes. Even though there is still a tendency that the baryon asymmetry is erased by 
the chemical equilibrium due to the sphaleron transitions, the finite mass of the tau lepton 
shifts the chemical equilibrium between B and L 3 towards the B side and leaves a finite 
asymmetry in the end. Their estimate is 




B = -^>: [Ih--B ! + _-£ (207) 



where the temperature T ~ Tc ~ 200 GeV is when the sphaleron transition freezes out 
(similar to the temperature of the electroweak phase transition) and m T (T) is expected to 
be somewhat smaller than m T (0) = 1.777 GeV. Overall, the sphaleron transition suppresses 
the baryon asymmetry by a factor of ~ 10~ 6 . This suppression factor is sufficient to keep the 
total baryon asymmetry at a reasonable order of magnitude in many of the cases discussed 
above. 



6 Dark Matter and Accelerator Constraints 

There is considerable evidence for dark matter in the Universe [H2|. The best observational 



evidence is found on the scale of galactic halos and comes from the observed flat rotation 
curves of galaxies. There is also good evidence for dark matter in elliptical galaxies, as 
well as clusters of galaxies coming from X-ray observations of these objects. In theory, we 
expect dark matter because 1) inflation predicts Q — 1, and the upper limit on the baryon 
(visible) density of the Universe from big bang nucleosynthesis is Qb < 0.1 [j§3]] ; 2) Even 
in the absence of inflation (which does not distinguish between matter and a cosmological 
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constant), the large scale power spectrum is consistent with a cosmological matter density of 
Q ~ 0.3, still far above the limit from nucleosynthesis; and 3) our current understanding of 
galaxy formation is inconsistent with observations if the Universe is dominated by baryons. 

It is also evident that not only must there be dark matter, the bulk of the dark matter 
must be non-baryonic. In addition to the problems with baryonic dark matter associated with 
nucleosynthesis or the growth of density perturbations, it is very difficult to hide baryons. 
There are indeed very good constraints on the possible forms for baryonic dark matter in 
our galaxy. Strong cases can be made against hot gas, dust, jupiter size objects, and stellar 



remnants such as white dwarfs and neutron stars [g4[ . 

In what follows, I will focus on the region of the parameter space in which the relic 
abundance of dark matter contributes a significant though not excessive amount to the 
overall energy density. Denoting by Q x the fraction of the critical energy density provided 
by the dark matter, the density of interest falls in the range 

0.1 < Q x h 2 < 0.3 (208) 

2 



The lower limit in eq. (|208| ) is motivated by astrophysical relevance. For lower values of Q x h 



there is not enough dark matter to play a significant role in structure formation, or constitute 



a large fraction of the critical density. The upper bound in ( |208|) , on the other hand, is an 



absolute constraint, derivable from the age of the Universe, which can be expressed as 

H t = J dx (l - ft - Q A + Q A x 2 + fi/x) ~ V2 (209) 



In ( |209|) , Q is the density of matter relative to critical density, while fl A is the equivalent 



contribution due a cosmological constant. Given a lower bound on the age of the Universe, 
one can establish an upper bound on Qh 2 from eq.( [209| ). A safe lower bound to the age of 



the Universe is t > 12 Gyr, which translates into the upper bound given in ( |208| ). Adding 
a cosmological constant does not relax the upper bound on Qh 2 , so long as Q + Q\ < 1. If 
indeed, the indications for a cosmological constant from recent supernovae observations ME 



turn out to be correct, the density of dark matter will be constrained to the lower end of the 
range in (|20£ 



As these lectures are focused on supersymmetry, I will not dwell on the detailed evidence 
for dark matter, nor other potential dark matter candidates. Instead, I will focus on the role 
of supersymmetry and possible supersymmetric candidates. As was discussed at the end of 
section 3, one way to insure the absence of unwanted, B and L- violating superpotential terms, 
is to impose the conservation of .R-parity. In doing so, we have the prediction that the lightest 
supersymmetric particle (LSP) will be stable. It is worth noting that -R-parity conservation 
is consistent with certain mechanisms for generating neutrino masses in supersymmetric 
models. For example, by adding a new chiral multiplet u c , along with superpotential terms 
of the form, H2Lu c + Mz/V c , although lepton number is violated (by two units), i?-parity is 
conserved. In this way a standard see-saw mechanism for neutrino masses can be recovered. 

The stability of the LSP almost certainly renders it a neutral weakly interacting particle 
I9| . Strong and electromagnetically interacting LSPs would become bound with normal 



matter forming anomalously heavy isotopes. Indeed, there are very strong upper limits on 
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the abundances, relative of hydrogen, of nuclear isotopes fSSfl , n/nn & 10 15 to 10 29 for 



1 GeV < m < 1 TeV. A strongly interacting stable relics is expected to have an abundance 
n/njj ^ 10~ 10 with a higher abundance for charged particles. 

There are relatively few supersymmetric candidates which are not colored and are elec- 
trically neutral. The sneutrino |87| is one possibility, but in the MSSM, it has been excluded 
as a dark matter candidate by direct |88| searches, indirect |89] and accelerator |H) searches. 
In fact, one can set an accelerator based limit on the sneutrino mass from neutrino counting, 
mo ?t 43 GeV PI| . In this case, the direct relic searches in underground low-background 
experiments require > 1 TeV [9^] . Another possibility is the gravitino which is probably 



the most difficult to exclude. I will concentrate on the remaining possibility in the MSSM, 
namely the neutralinos. 

The neutralino mass matrix was discussed in section 3.3 along with some particular 
neutralino states. In general, neutralinos can be expressed as a linear combination 



X = aB + f3W 3 + 7 #! + SH 2 



(210) 



and the coefficients a,/3, 7, and 5 depend only on M 2 , //, and tan/3 (assuming gaugino mass 
unification at the GUT scale so that Mi = f — Mo). 



There are some limiting cases in which the LSP is nearly a pure state | I9|| . When // — > 0, 
S° is the LSP with 

^^/z = /isin2/3 (211) 
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When M 2 -> 0, the photino is the LSP with E3| 
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When M 2 is large and M 2 <C \i then the bino B is the LSP [^3| and 



77lg ~ 



Mi 
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jj, for /i > are the LSPs depending on the sign of ji 



and finally when fi is large and fi M 2 the Higgsino states i?(i2) with mass m H {12) — ~^ 
for n < 0, or if[ 12 ] with mass 

In Figure^ [p3| , regions in the M 2 ,/i plane with tan/3 = 2 are shown in which the LSP 
is one of several nearly pure states, the photino, 7, the bino, B, a symmetric combination 
of the Higgsinos, H^ 2 ), or the Higgsino, S. The dashed lines show the LSP mass contours. 
The cross hatched regions correspond to parameters giving a chargino (W ± ,H ± ) state with 



mass rriy < A5GeV and as such are excluded by LEPK 



by LEP1.5 



and LEP2 |6 



This constraint has been extended 
and is shown by the light shaded region and corresponds 



to regions where the chargino mass is < 95 GeV. The newer limit does not extend deep 
into the Higgsino region because of the degeneracy between the chargino and neutralino. 
Notice that the parameter space is dominated by the B or H\ 2 pure states and that the 
photino (often quoted as the LSP in the past |23|, |97|) only occupies a small fraction of the 
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Figure 9: Mass contours and composition of nearly pure LSP states in the MSSM [93 . 

parameter space, as does the Higgsino combination S°. Both of these light states are now 
experimentally excluded. 

The relic abundance of LSP's is determined by solving the Boltzmann equation for the 
LSP number density in an expanding Universe. The technique ||98|| used is similar to that 
for computing the relic abundance of massive neutrinos p9|. The relic density depends on 
additional parameters in the MSSM beyond M 2 ,/i, and tan/3. These include the sfermion 
masses, mj, the Higgs pseudo-scalar mass, m^, and the tri-linear masses A as well as two 
phases 6^ and 6 a- To determine, the relic density it is necessary to obtain the general 
annihilation cross-section for neutralinos. This has been done in [ |100| , [TUT], |102| , |103[ |. In 
much of the parameter space of interest, the LSP is a bino and the annihilation proceeds 
mainly through sfermion exchange as shown in Figure 10. For binos, as was the case for 



photinos 1^3], |9?|, it is possible to adjust the sfermion masses to obtain closure density in a 
wide mass range. Adjusting the sfermion mixing parameters [|104| or CP violating phases 
34|, |35| allows even greater freedom. 





f 

Figure 10: Typical annihilation diagram for neutralinos through sfermion exchange. 
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Because of the p-wave suppression associated with Majorana fermions, the s-wave part 
of the annihilation cross-section is suppressed by the outgoing fermion masses. This means 
that it is necessary to expand the cross-section to include p-wave corrections which can be 
expressed as a term proportional to the temperature if neutralinos are in equilibrium. Unless 
the B mass happens to lie near mz/2 or in which case there are large contributions to 

the annihilation through direct s-channel resonance exchange, the dominant contribution to 
the BB annihilation cross section comes from crossed t-channel sfermion exchange. In the 
absence of such a resonance, the thermally-averaged cross section for BB —>■ ff takes the 
generic form 



(av) = {l-^r) 1 ' 2 ^- 
v 1 v m\> 128vr 



= a + bx (214) 

where Y^pt) is the hypercharge of Jl{r), A/ = m 2 ~ + m| — mj, and we have shown only the 
leading P-wave contribution proportional to x = T/m^. Annihilations in the early Universe 
continue until the annihilation rate T ~ avn x drops below the expansion rate, H. For 
particles which annihilate through approximate weak scale interactions, this occurs when 
T ~ m x /20. Subsequently, the relic density of neutralinos is fixed relative to the number of 
relativistic particles. 

As noted above, the number density of neutralinos is tracked by a Boltzmann-like equa- 
tion, 

d i t =-^n-(av){n*-nl) (215) 

where no is the equilibrium number density of neutralinos. By defining the quantity / = 
n/T 3 , we can rewrite this equation in terms of 

df / 1 \ 1/2 



dx m A^ KN ) v'-® (216) 

The solution to this equation at late times (small x) yields a constant value of /, so that 
n oc T 3 . The final relic density expressed as a fraction of the critical energy density can be 



written as [|19|] 




a x h 2 ~ 1.9 x 10- n ^ Ni /A ' 2 (217) 



where (T x /T 7 ) 3 accounts for the subsequent reheating of the photon temperature with respect 
to Xi due to the annihilations of particles with mass m < Xfin x [j4^| . The subscript / refers 
to values at freeze-out, i.e., when annihilations cease. 

In Figure |ll] ||105|| , regions in the M 2 — /i plane (rotated with respect to Figure [J) with 



tan/3 = 2, and with a relic abundance 0.1 < Qh 2 < 0.3 are shaded. In Figure O, the 



sfermion masses have been fixed such that m = 100 GeV (the dashed curves border the 
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region when mo = 1000 GeV). Clearly the MSSM offers sufficient room to solve the dark 
matter problem. In the higgsino sector if 12 , additional types of annihilation processes known 
as co- annihilations [106, 108 , 107|| between #(12) an d the next lightest neutralino (-ff[i2]) must 
be included. These tend to significantly lower the relic abundance in much of this sector and 
as one can see there is little room left for Higgsino dark matter ||105| . 

As should be clear from Figures ^| and 11, binos are a good and likely choice for dark 
matter in the MSSM. For fixed mj, Qh 2 > 0.1, for all = 20 — 250 GeV largely independent 
of tan/? and the sign of /i. In addition, the requirement that m^ > translates into an 
upper bound of about 250 GeV on the bino mass |93| . |109|| . By further adjusting the trilinear 
A and accounting for sfermion mixing this upper bound can be relaxed ||104|| and by allowing 
for non-zero phases in the MSSM, the upper limit can be extended to about 600 GeV [|34 |. 
For fixed Qh 2 = 1/4, we would require sfermion masses of order 120 - 250 GeV for binos 
with masses in the range 20 - 250 GeV. The Higgsino relic density, on the other hand, is 
largely independent of mj. For large /i, annihilations into W and Z pairs dominate, while 
for lower /i, it is the annihilations via Higgs scalars which dominate. Aside from a narrow 
region with m^ i2 < m\y and very massive Higgsinos with m^ i2 > 500 GeV, the relic density 
of H12 is very low. Above about 1 TeV, these Higgsinos are also excluded. 




Figure 11: Regions in the M 2 -/x plane where 0.1 < Qh 2 < 0.3 [105 . Also shown are the 
Higgsino purity contours (labeled 0.1, 0.5, and 0.9). As one can see, the shaded region is 
mostly gaugino (low Higgsino purity). Masses are in GeV. 

As discussed in section 4, one can make a further reduction in the number of parameters 
by setting all of the soft scalar masses equal at the GUT scale, thereby considering the 
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CMSSM. For a given value of tan (3, the parameter space is best described in terms of the 



common gaugino masses and scalar masses tuq at the GUT scale. In Figure [T^ [35fl, this 
parameter space is shown for tan/3 = 2. The light shaded region corresponds to the portion 
of parameter space where the relic density Q x h 2 is between 0.1 and 0.3. The darker shaded 
region corresponds to the parameters where the LSP is not a neutralino but rather a tr. In 
the rriQ — m^j plane, the upper limit to m is due to the upper limit Q x h 2 < 0.3. For larger 
mo, the large sfermion masses suppress the annihilation cross-section which is then too small 
to bring the relic density down to acceptable levels. In this region, the LSP is mostly B, 
and the value of fi can not be adjusted to make the LSP a Higgsino which would allow an 
enhanced annihilation particularly at large itlq. The cosmo logically interesting region at the 
left of the figure is rather jagged, due to the appearance of pole effects. There, the LSP can 
annihilate through s-channel Z and h (the light Higgs) exchange, thereby allowing a very 
large value of m . Because the fL./t 2 = 0.3 contour runs into the f^-LSP region at a given 



value of mi/2, it was thought [[32] that this point corresponded to an upper limit to the LSP 
mass (since is approximately 0.4 m^. As we will see, this limit has been extended due 
to co-annihilations of the B and tr [ 110|| . 




Figure 12: Region in the mi/2-mo plane where 0.1 < Qh 2 < 0.3 P5] . Masses are in GeV. 

The mi/2 — mo parameter space is further constrained by the recent runs at LEP. The 
negative results of LEP 1 searches for Z° — > x + X~ an d Z° — » xx' (where x' denotes a 
generic heavier neutralino) already established important limits on supersymmetric model 



parameters, but left open the possibility that the lightest neutralino might be massless [94 



Subsequently, the data from higher-energy LEP runs, based on chargino and neutralino pair 
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production, complemented LEP 1 data and excluded the possibility of a massless neutralino, 
at least if the input gaugino masses M a were assumed to be universal ||111| , ^1 



In Figure 13 [112 



slepton searches [111 



the constraints imposed by the LEP chargino and neutralino and 
jf| (hatched region) at LEP 2 are shown. The hatched regions corre- 
spond to the limits in the MSSM. The thick curve to the right of this region, corresponds to 
the slightly more restrictive bound due to the assumption of universal Higgs masses at the 
GUT scale (CMSSM). The distinction between these two cases is more apparent at other 
values of tan/3 pT] , |112j| . The bounds shown here correspond to the run at LEP at 172 GeV 
center of mass energy. The DO limit on the gluino mass |113|| is also shown (dotted line). 
Of particular importance are the bounds on supersymmetric Higgs production, for which 
we consider both the e + e~ — > hZ and e + e~ — > hA reactions. The regions bounded by the 
lack of Higgs events are labeled nUHM corresponding to the MSSM (non-Universal Higgs 
Mass) and UHM corresponding to the UHM. Once again, the limits plotted in this figure 
(|TB|) correspond to the 172 GeV run and have been greatly improved since. In this case the 
improvement in the bound when restricting the value of fi to take its CMSSM value is clear. 
At lower tan/3, the constraint curve moves quickly to the left, i.e., to higher values of 




Figure 13: The domains of the ^1/2^^0 plane (masses in GeV) that are excluded by the 
LEP 2 chargino and selectron searches, both without (hatched) and with the assumption of 
Higgs scalar-mass universality for fi < and tan /3 — 2 [ |112|| . Also displayed are the domains 
that are excluded by Higgs searches (solid lines) with and without the assumption of universal 
scalar masses for Higgs bosons (UHM). The region marked theory is excluded cosmologically 
because viif R < m x . The domains that have relic neutralino densities in the favored range 
with (dark) and without (light shaded) the scalar-mass universality assumption. 

The region of the mi/2 — itlo plane in which 0.1 < Q x h 2 < 0.3 for some experimentally 
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allowed value of /i < is light-shaded in Figure [13], and the region of the plane in which 
0.1 < £l x h 2 < 0.3 for fi determined by the CMSSM constraint on the scalar masses is shown 
dark-shaded. The MSSM region extends to large values of itlq since /i can be adjusted to 
take low values so the the LSP is Higgsino like and the relic density becomes insensitive to 
the sfermion masses. 

As one can see from Figure [13], the combined bounds from the chargino and slepton 
searches, provide a lower bound to m\/2 which can be translated into a bound on the neu- 
tralino mass. In the region of interest, we have the approximate relation that m x « 0.4t%/2. 



200 



ioo- 



m _j=95 

r 



a) tan p= 3, \i < 



200 




1r % 



100 200 300 400 500 600 
m i/2 




21)0 



Si c)tanp=10, h<0 




200- 



100- 



d) tanp=10, H>0 



m =95 
X 



100 200 300 400 500 600 

m i/2 



100 200 300 400 500 600 



Figure 14: The light-shaded area is the cosmologically preferred region with 0.1 < Qh 2 < 0.3. 
The light dashed lines show the location of the cosmologically preferred region if one ignores 
coannihilations with the light sleptons. In the dark shaded regions in the bottom right of 
each panel, the LSP is the tr, leading to an unacceptable abundance of charged dark matter. 
Also shown are the isomass contours m y ± = 95 GeV and = 95, 100, 105, 110 GeV, as well 



as an indication of the slepton bound from LEP [ |115|| . These figures are adapted from those 



in |TT0 |. 



Ultimately, the Higgs mass bound alone does not provide an independent bound on 
mi/2- The reason is that the Higgs constraint curves bend to the left at large mo, where 
large sfermion masses lead to greater positive radiative corrections to the Higgs mass, and 
the Higgs curve strikes the chargino bound at sufficiently large itlq. However, when combined 
with cosmological limit on the relic density, a stronger constraint can be found. Recall that 
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at large m , the cosmological bound is satisfied by lowering \i. At certain values of tan/3 
and rrii/2, one can not lower fi and remain consistent with the relic density limit and the 
Higgs mass limit. This is seen in Figure ^ by the short horizontal extension of the nUHM 
Higgs curve. At lower values of tan/3 this extension is lengthened. In the UHM case, once 
again the UHM Higgs curve bends to the right at large mo. However, cosmology excludes 
the large m region in the CMSSM (CMSSM and UHM are treated synonymously here). In 
this case the mass bound on m 1( / 2 or m x is much stronger. 

The results from the 172 GeV run at LEP provided a bound on the chargino mass of 
m x ± > 86 GeV and corresponding bound on the neutralino mass of m x > 40 GeV [|112 |. 
As noted above, the UHM Higgs mass bound becomes very strong at low tan/3. In fact, 
for tan/3 < 1.7, the UHM Higgs curve moves so far to the left so as to exclude the entire 
dark shaded region required by cosmology. Subsequent to the 183 GeV run at LEP ||114|| , 



the chargino mass bound was pushed to m x ± > 91 GeV, and a Higgs mass bound was 
established to be m h > 86 GeV for tan /3 < 3, and m h > 76 GeV for tan /3 > 3. These limits 
were further improved by the 189 GeV run, so that at the kinematic limit the chargino mass 
bound is m x ± > 95 GeV, implying that m 1( / 2 > 110 GeV and the neutralino mass limit 
becomes m x > 50 Gev. The Higgs mass bound is now (as of the 189 run) m h > 95 GeV for 



low tan (3 115 
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Figure 15: The same as Fig. Ill, for /i < 0, but extended to larger mi/ 2 . These figures are 



adapted from those in [110 . 



In the MSSM, in the region of the M 2 -fi plane, where -^(12) is the LSP (that is, at large fi 
and very large M 2 ), the next lightest neutralino or NLSP, is the H^ 2 ] state and the two are 
nearly degenerate with masses close to fi. In this case, additional annihilation channels (or 
co-annihilation) which involve both a -^(12) and a -^[12] in the initial state become important 
106| , |107|| . The enhanced annihilation of Higgsinos lowers the relic density substantially 



and virtually eliminates the Higgsino as a viable dark matter candidate. In the CMSSM, 
co-annihilation is also important, now at large values of m 1( / 2 [ |110]| . Recall, that previously 
we discussed a possible upper limit to the mass of the neutralino in the CMSSM, where the 



cosmologically allowed region of Figure [12] runs into the region where the tr is the LSP. Along 
the boundary of this region, the neutralino B, and the tr are degenerate. Therefore close to 
this boundary, new co-annihilation channels involving both B and tr (as well as the other 
right-handed sleptons which are also close in mass) in the initial state become important. 
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Figure 16: As in [bj for tan (3 = 2, 2.5, and 3, for both \i > and fi < 0. Higgs mass contours 
of 86, 95, and 96 GeV, are displayed to show the dependence on tan/5. I thank Toby Falk 
for providing these figures. 
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As in the MSSM, the co-annihilations reduce the relic density and as can be seen in Figures 
14| the upper limit to the neutralino is greatly relaxed ||1 10|| . As one can plainly see, the 



cosmologically allowed region is bent way from the f^-LSP region extending to higher values 
of mi/2- Also shown in this set of figures are the iso-mass contours for charginos, sleptons 
and Higgses, so that the LEP limits discussed above can be applied. In these figures, one 
can see the sensitivity of the Higgs mass with tan (3. 

Despite the importance of the coannihilation, there is still an upper limit to the neutralino 
mass. If we look at an extended version of Figures H, as shown in Figures 15, we see that 
eventually, the cosmologically allowed region intersects the f^-LSP region. The new upper 
limit occurs at m 1( / 2 ~ 1450 GeV, implying that m x < 600 GeV. 

As noted above, the Higgs mass bounds can be used to exclude low values of tan/3. 
However, when co-annihilations are included, the limits on tan/3 is weakened. In the m^j — 
tuq plane, the Higgs iso-mass contour appears nearly vertical for low values of m . A given 
contour is highly dependent on tan (3, and moves to the right quickly as tan (3 decreases. 
This behavior is demonstrated in the series of Figures [H], which show the positions of the 
Higgs mass contours for tan/? = 2,2.5, and 3, for both positive and negative /i. If we 
concentrate, on the 95 GeV contour, we see that at tan (3 = 3, the bulk of the cosmological 
region is allowed. At tan/3 = 2.5, much of the bulk is excluded, though the trunck region is 
allowed. At the lower value of tan (3 = 2, the 95 GeV contour is off the scale of this figure. A 
thorough examination |110|1 yields a limit tan/3 > 2.2 for /i < 0. For positive /i, the results 
are qualitatively similar. The Higgs mass contours are farther to the left (relative to the 
negative /i case), and the limit on tan/3 is weaker. Nevertheless the limit is tan/3 > 1.8 for 
/i > 0. 

As the runs at LEP are winding down, the prospects for improving the current limits or 
more importantly discovering supersymmetry are diminishing. Further progress will occur in 
the future with Run II at the Tevetron at ultimately and the LHC. Currently, while we have 
strong and interesting limits on the MSSM and CMSSM parameter space for LEP, much of 
the phenomenological and cosmological parameter space remains viable. 
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